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ABSTRACT 

Low-resolution spectra from 3000-9000 A of young low-mass stars and brown dwarfs were obtained 
with LRIS on Keck I. The excess UV and optical emission arising in the Balmer and Paschen continua 
yields mass accretion rates ranging from 2 x 10~^^ to 10~^ Mq yr~^. These results are compared with 
HST/STl'^ spectra of roughly solar-mass accretors with accretion rates that range from 2 x 10"^" to 
5 X 10^^ M© yr^^. The weak photospheric emission from M-dwarfs at < 4000 A leads to a higher 
contrast between the accretion and photospheric emission relative to higher-mass counterparts. The 
mass accretion rates measured here are systematically ~ 4 — 7 times larger than those from Ha 
emission line profiles, with a difference that is consistent with but unlikely to be explained by the 
uncertainty in both methods. The accretion luminosity correlates well with many line luminosities, 
including high Balmer and many He I lines. Correlations of the accretion rate with Ha 10% width 
and line fluxes show a large amount of scatter. Our results and previous accretion rate measurements 
suggest that M oc M^-^''^^-'^^ for accretors in the Taurus Molecular Cloud. 
Subject headings: stars: pre-main sequence — stars: planetary systems: protoplanetary disks 



1. INTRODUCTION 

In the magnetospheric accretion model for accretion 
onto young stars, the stellar magnetosphere truncates 
the disk at a few stellar radii. Gas from the disk accretes 
onto the star along magnetic field lines and shocks at 
the stellar surface. The ^ 10* K optically-thick post- 
shock and optically-thin pre-shock gas produce emission 
in the Balmer and Paschen continuum and in many lines, 
including the B almer and Paschen series and the Ca II IR 
triplet (Hartmann et al.lll994[Calvet fc Gullbring|ll998f ). 
The accretion luminosity may be estimated from any of 
these accretion tracers and subsequently converted into 
a mass accretion rate. 

Historically, mass accretion rates onto young stars 
have been measured directly using optical veiling from 
the Paschen continuum and stronger Balmer contin- 
uum emission shortward of the Balmer jump. Mod- 
elling of the Paschen and Balmer continuum emission 
with a generic isothermal, plane-parallel slab yields M ~ 
10""^ — 10~^ ° Mf7) yr~^ for accretors with masses o f 
0.4-1.0 M© (iValenti et al.l fl99l iGu llbring et al.l ll99l. 
iCalvet fc G uUbrind (|1998h obtained similar accretion 
rates by applying more realistic shock geometries with 
a range of temperatures to explain the observed ac- 
cretion continuum. The accretion luminosity can be 
measured instead at longer wavelengths with opti- 
cal veiling measurements from high-dispersion spectra 
(e.g. [Ba sri & Batallia 1990; Hartigan et al. 1991, 1995; 
iHartigan fc KenvonI 120031 : IWhite fc 'Hillenbrandl l2004h . 
albeit with a larger uncertainty in bolometric correction. 

Several deep surveys have since identified large popu- 
lations of young very low-mass stars and brown dwarfs 
with inasses < 0.2 M© (e.g.lB riccno et al. 1998; Luhman 
[20041: iSlesnick et al.l[200lTGuieu et al.ll2007n . Excess IR 
emission indicates the presence of disks around some of 
these very low-mass objects^ (e.g. iMuzerolle et al.l[2003t 
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iJavawardhana et all 12001 iNatta et ani2004[ ). The ex- 
cess Balmer continuum emission produced by accretion 
onto very low-mass stars and brown dwarfs has not been 
observed previously. Instead, accretion has been been 
measured based on magnetospheric models of the Ha 
hue profile (Muzerolle et al. 2000, 2003, 2005), by 
extending relationships between Al and line luminosi- 
ties established at higher masses to brown dwarf masses 
(Natta et al. 2004, 2006; Mohanty et al. 2005), or in a 
few cases from measu rements of the Paschen continuum 
(jWhite fc Basril[2003h . 

In this paper, we present the first observations of the 
Balmer jump, and hence direct measurements of the ac- 
cretion luminosity, for young very low-mass stars and 
brown dwarfs using data obtained with LRIS on Keck I. 
We compare the Balmer and Paschen continuum emis- 
sion from very low-mass stars to that from higher-mass 
accretors that were observed with HST. The weak U- 
band photospheric emission of M-dwarfs leads to a high 
contrast between accretion and photospheric emission at 
A < 3700 A. The properties of the Balmer continuum 
and jump for very low-mass stars and brown dwarfs are 
mostly similar to those of the higher-mass accretors, al- 
though the sensitivity to low M improves to smaller stel- 
lar mass. The Balmer jump tends to be larger for smaller 
mass accretion rates, indicating lower optical depth for 
the emitting gas. The luminosity in many lines, includ- 
ing the Balmer series and He I lines, correlate well with 
the accretion luminosity. Our A'l estimates are system- 
atically ~ 4 — 7 times higher than those estimated from 
modelling the Ha line profile. We find that M oc Af^-^ 
for accretors in Taurus and discuss biases in the selec- 
tion of those accretors that may affect this calculated 
relationship. 

2. OBSERVATIONS 

^ The objects studied in this paper span the mass range from 
0.024—0.2 Mq\ hereafter all are referred to as "stars" for ease of 
terminology. We refer to ^ 0.3 — 1 Mq stars as "higher-mass stars." 
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2.1. Keck I/LRIS 

We used LRIS (|Oke et al.lll995l : lMcCarthv et al.lll998D 

on Keck I to obtain low-resolution optical spectra cover- 
ing 3000-9000 A of 17 young, low-mass Taurus sources 
and 2M1207-3932 on 23 November 2006. We observed 
2M0436-f 2351 and obtained a deeper observation of 
2M1207-3932 on 7 February 2007. Our observation log 
(Table 1) lists the exposure time for each target and ab- 
breviations for each target used throughout this paper. 

Most targets were observed with a 175" x I'.'O slit. The 
light was split into separate red and blue beams with the 
D560 dichroic. The red light illuminates a 2048 x 2048 
CCD with 0'.'21 pixels and the blue light illuminates 
a 2048 X 4096 CCD with 0'.'135 pixels. We used the 
400/8500 grism and the GG495 filter in the red chan- 
nel and the 400/3400 grism in the blue channel. Each 
red pixel covers ^ 1.85 A, yielding a resolution R ~ 1400 
over the 5400-9000 A wavelength range. Each blue pixel 
covers ~ 1.1 A, yielding a resolution R ^ 700 over the 
3000-5700 A wavelength range. For each target we typ- 
ically obtained multiple red spectra during a 1-2 long 
blue exposure. 

We bias-subtracted the 23 November blue exposures 
and estimated a bias subtraction for the 7 February 2007 
blue exposure from images with no blue emission. An 
overscan subtraction was applied to all red images. The 
red images were flat-fielded with dome flats, the blue 
images from November 2006 were flat-flelded with sky 
flats, and the blue images from February 2007 were flat- 
fielded with halogen lamp flats. We located each spec- 
trum on the detector as a function of position in the dis- 
persion direction. The counts in the red and blue spectra 
were extracted with 11 and 13-pixel widths, respectively. 
The sky background was measured from nearby windows 
and subtracted from the spectrum. The spectral regions 
within 25 A of 4590 and 5270 A include a transient ghost 
and are ignored. 

The relative wavelength solution was calculated from a 
HgNeArCdZn lamp spectrum. The absolute wavelength 
solution for each observation was calculated by measur- 
ing the positions of sky emission lines. 

All stars were observed at the parallactic angle. On 23 
November the seeing improved from 0'.'9 to 0'.'7 from the 
beginning to the middle of the night and remained stable 
at 0'.'7 for the second half of the night. On 7 February the 
seeing was 0'.'7 during our observation of 2M0436-I-2351 
and 1'.'2 during our observation of 2M1207-3932. 

We obtained observations of the spectropho t omet- 
ric standards (iMasse y et aH 119881 : lOke et all Il990( ) 
PG0205-hl34 and PG0939+262 at the beginning and 
end of the November 2006 run and G191B2B and 
PG1545-h035 on February 7 2007 with the I'.'O aperture. 
The November 23 standards have published fluxes from 
3100-8000 A. We approximate the flux of PG0939-I-262 
from 8000-9000 A by assuming that the white dwarf 
spectrum is a blackbody, and confirm this calibration 
by comparing the spectra of 2M1207-3932 obtained on 
both nights. The atmospheric extinction correction was 
obtained from multiple observations of G191B2B on 7 
February 2007. 

Slit losses were ~ 70% at the beginning of the night and 
~ 22% at the end of the night in Nov. 2006. In adjacent 
red spectra of the same target, the extinction-corrected 



count rate increased early in the night and remained sta- 
ble during the second half of the night, mirroring similar 
changes in the seeing. This increase in count rate can 
account entirely for the factor of ~ 3 in slit losses. The 
count rates in the first half of the night are therefore 
multiplied by a factor between 1 — 3, based on when the 
observation occurred. The small wavelength dependence 
in this factor is found by interpolating from the observa- 
tions of the photometric standards. 

Table 1 lists synthetic UBVR magnitudes for each 
source and our estimated uncertainty in absolute flux 
calibration. The estimated photometry is within 0.1-0.5 
mag of most existing literature measurements, with the 
exception of several strong accretors. The initial flux 
calibration of 2M1207-3932 on the two nights differed by 
0.27 mag. We estimate that the absolute fluxes are cali- 
brated to at worst a factor of 2 early in the night and a 
factor of 1.3 in the second half of the night for the Nov. 
2006 observations. 

2.2. Archival HST Spectra 

We supplement our new Keck data with existing blue 
775T/STIS spectra of late-K or early-M accretors and 
two young non-accretors taken from the MAST archive 
(see Table 2). The spectra were obtained with the G430L 
grating and span 2900-5700 A with a spectral resolution 
of ^ 600. Several of these observations were discussed in 
Herczeg et al. (2004, 2005, 2006), including calculations 
of the mass accretion rate. We re-analyze these data 
here to provide a small set of higher-mass accretors with 
spectra analyzed in a manner consistent with our Keck 
//LRIS spectra to compare with results from our sample 
of low-mass accretors. 

2.3. Description of our sample 

Several targets (CM Tau, CIDA 1, 2M0414-f 2811) in 
our Keck //LRIS program were selected based on ex- 
pectation of a high mass accretion rate inferred from 
strong excess continuum emission and optical line emis- 
sion. Other targets were selected based on a combina- 
tion of spectral type, low extinction, existing M mea- 
surements from Ha, and visual magnitude. MHO 5 and 
KPNO 11 were selected also for a full- width of the Ha 
line at 10% the p eak flux (h ereafter Ha 10% wid th) of 
150-2 00 km s'^ (|Muzerolle et al. 2003; Mohant v et al.l 
l2005f ). which is intermediate between definite accretors 
and definite non-accretors for their respective spectral 
types. Three relatively-bright young stars in Taurus were 
chosen as photospheric templates based on spectral type 
and lack of strong accretion, as identified by by weak 
Ha emission and absence of IR excess emission. Since 
the template stars are young, their surface gravity will 
be similar to the rest of the sample and provide a bet- 
ter match to the photospheric spectrum than could be 
obtained from older field dwarfs. 

To obtain a larger range in M for analyzing the prop- 
erties of the accretion continuum and relationships be- 
tween line and continuum emission (see §5 and 6), we 
supplemented the Keck //LRIS spectra with //ST/STIS 
observations of stars with higher central masses. This 
sample consists of most available HST observations of 
the 3000-5000 A region for accreting young stars with 
spectral types late K or early M. 
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Table 3 describes the properties of these targets (see 
also §9). Muhiple LRIS observations of 2M1207-3932 
and HST observations of DF Tau and TW Hya are 
treated independently throughout this paper. 

3. OBSERVATIONS OF ACCRETION DIAGNOSTICS 

Figure 1 compares three M5-M5.5 stars in our 
Keck/l.mS> spectra fro m 3200-5600 A. Strong T iO and 
CaH molecular bands (|Kirkpatrick at ahl Il993f ) domi- 
nate the photospheric emission of mid-M dwarfs. At 
> 3700 A, such features are strong from V927 Tau 
and 804394-2336 but shallow from S05184-2327. This 
veiling of the photospheric emission from S0518-I-2327 
is attributed to excess Paschen continuum emission. 
S0439-h2336 and S0518-h2327 both exhibit strong ex- 
cess emission shortward of jump at 3700 A, which is 
attributed to excess Balmer continuum emission. Line 
emission, including in the Balmer series, is much stronger 
from S0439-h2336 and 80518-^2327 than from V927 Tau. 
The excess Balmer and Paschen continuum and emission 
lines detected from young stars are commonly attributed 
to accretion. In the following subsections we describe 
and measure these accretion diagnostics for each source. 

3.1. Measuring the Paschen Continuum 

Veiling"^ is measured by comparing the depth of ab- 
sorption features, which are filled in by emission from the 
accretion continuum, to photospheric templates. These 
features include the TiO and CaH absorption bands and 
also the strong Ca I A4226 absorption line. The Ca I 
line depth may be related to chromos pheric structure 
and may be affected by stellar rotation (jMauas fc Falchil 
[1996; Short ct al. 1997). 

The young M-dwarfs V927 Tau (M5), MHO 7 (M5.5), 
CFHT 7 (M6.5) and 2M0439-f 2544 (M7.25) serve as pho- 
tospheric templates (see §9.1). These stars provide a 
better match to the surface gravity and chromospheric 
emission of young accretors than would older field M- 
dwarfs. 2M0439-I-2544 is not an ideal template because 
of ongoing accretion, and V927 Tau may also be accret- 
ing (see §3.4 and §9.4.2). Use of 2M0439-h2544 and V927 
Tau as templates may cause us to underestimate the veil- 
ing by 0.05 and 0.02, respectively, at 5000 A and 0.01 
each at 7000 A. At < 4500 A we use CFHT 7 instead 
of 2M0439+2544 as a photospheric template for late M- 
stars. 

Large veilings (> 2) are difficult to measure because 
the high ratio of accretion-to-photosphcric flux masks 
any photospheric absorption features. However, a large 
error for a high value of veiling leads to only a small 
error in the measurement of the underlying continuum 
flux. Our lower detection limit of veiling measurements 
is ^ 0.05 — 0.1, depending on how well the photospheric 
template matches the spectrum. The veiling is calculated 
from regions that do not include obvious strong emission 
lines. For weaker accretors, any such lines may fill in 
some of the absorption and thereby artificially increase 
our veiling estimate for some stars, particularly at 4227 
A (see Fig. 2). 

Figures 2-4 show three spectral regions where veiling 
is measured. Table 4 lists veiling measurements and the 

^ The veiling at wavelength A, r\, is defined by the accretion 
continuum flux {Fobaerved - Fphotosphere) divided by the photo- 
spheric flux. 



corresponding fiux in the accretion continuum at dis- 
tinct wavelengths between 4000-9000 A. Excess contin- 
uum emission is detected longward of 4000 A for many 
stars in our sample and is identified as the Paschen con- 
tinuum. The accretion continuum is relatively fiat at 
< 6000 A for most targets, and increases at A > 6000 A 
for most sources where it is detected. The continuum flux 
from 2M0414-H2811 at 8600 A is substantially lower than 
that at 7800 A, which may be consistent with weaker 
Brackett continuum emission longward of the Paschen 
jump. CIDA 1 shows no such flux reduction, which may 
suggest an alternate veiling source. 

The early-M and late-K stars from the STIS data do 
not have the strong TiO and CaH absorption bands for 
veiling measurements. In several locations, the photo- 
spheric emission drops by a factor of 2-2.5 within ~ 200 
A. Table 5 fist veihngs at 4600 and 5400 A for this sam- 
ple, measured by comparing the depth of these features 
with photospheric templates, and at 4000 A, measured 
from the slope of the photospheric emission. We detect 
Paschen continuum emission from DF Tau, TW Hya, 
V836 Tau, DC Tau, and RU Lup. 

3.2. Measuring the Balmer Continuum 

When seen from young stars, excess Balmer continuum 
emission shortward of ~ 3700 A is typically attributed to 
accretion (Figure 5) . The Balmer limit occurs at 3646 A, 
but line blending in the Balmer series shifts the apparent 
jump to 3700 A. The observed Balmer jump, defined here 
as the ratio of fiux at 3600 A to the flux at 4000 A, ranges 
from ^ 0.35 for the non-accretors to 5.47 for the accretor 
2M1207-3932 (Table 6). 

Figure 6 compares spectra of stars with the weakest 
Balmer jumps in our sample to that of photospheric tem- 
plates from Table 4. Excess Balmer continuum emis- 
sion is clearly detected from 2M0455-I-3028, 80439-^2336, 
and V927 Tau relative to photospheric templates. Some 
weak, noisy Balmer continuum emission is also detected 
from MHO 5 and KPNO 11. MHO 7 and CFHT 7 may 
each show very mild increases in emission shortward of 
the Balmer jump, possibly from chromospheric emission 
(see §3.4). 

As defined here, the Balmer jump is the ratio of Balmer 
plus Paschen continuum flux at 3600 A divided by the 
Paschen continuum flux at 4000 A. The intrinsic Balmer 
Jump {J,nod in Table 6) is measured from the accretion 
continuum after dereddening the spectra (see §9.2) and 
subtracting a scaled photospheric spectrum. Table 6 also 
lists the measured slope of the continuum emission from 
3200-3650 A (Sb), scaled to the observed flux at 3600 

A. 

3.3. Emission Line Measurements 

Tables 7-10 list equivalent widths and fluxes for H 
Balmer lines. He lines, other strong lines associated with 
accretion, and lines usually associated with outflows.'* 
These tables are not complete linelists. We measure most 
emission line fluxes and equivalent widths by fitting a 
Gaussian profile to the observed line. The focus was 
poor at A > 8000 A and produced asymmetric emission 
lines. Line strengths in that spectral region are measured 
by summing the flux across a continuum-subtracted win- 
dow that includes the entire line. The Balmer and Ca II 
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H & K lines are bright from every source and are mea- 
sured directly from the spectra. For all other lines we 
subtract a scaled photospheric template to increase our 
sensitivity to small fluxes. These line equivalent widths 
and corresponding fluxes are relative to the photospheric 
template and may be offset by a small amount in equiva- 
lent width if line emission is present but undetected in the 
photospheric template. This method may also yield false 
weak detections if the depth of Na I or Ca II photospheric 
absorption features differ substantially from star-to-star. 
Weak Ca II IR triplet emission from several targets is 
only treated as real if all three of the lines are detected, 
unless the non-detection has a large upper limit on equiv- 
alent width. Our detection limits typically range from 
^ 0.1 — 1.5 A, depending on the match with a photo- 
spheric template and the signal-to-noise. The equivalent 
widths in low-resolution spectra are affected by absorp- 
tion and may differ from equivalent widths measured in 
high-resolution spectra because of different spectral re- 
gions used to calculate the continuum emission. 

The sources with strongest veiling also show the 
strongest line emission. For example, Ca II IR triplet 
emission is detected from only CIDA 1, 2M0414-t-2811, 
CM Tau, 2M0441-h2534, and S0518-f2327. Similarly, 
the Paschen series is detected from only CIDA 1, 
2M0414+2811, and CM Tau. Several Paschen lines blend 
with and contribute < 10% of the flux in the Ca II IR 
triplet lines. Other lines, including the Na I D doublet. 
He II A4686, and He I lines are detected from most stars 
with excess Balmer emission. 

Optical forbidden lines, including the [O I] A6300, 6363, 
[S II] A6716, 6731, and [S I] A4069, 4076 doublets, are de- 
tected from many sources. [N II] A6584 is detected from 
2M0444-H2612 but is difhcuh to resolve from the bright 
Ha emission for many other stars in the sample. The 
[O I] lines are unresolved in the cross-dispersion direc- 
tion, which constrains the emission to have a FWHM of 
< 0.2 - 0.3. 

3.4. Designation of Accretors from the Balmer 
continuum 

Many emission line diagnostics have been used to iden- 
tify accretion onto young very low-mass stars and brown 
dwarfs. An Ha equivalent width of > 10 A is the tra- 
ditional criterion above which a young star is classi- 
fied as an accretor. Against the weaker photospheric 
continua of redder M-dwarfs, chromospheric emission 
can produced Ha emission with equivalent widths larger 
than 10 A. In lieu of this diagnostic, the full-width at 
10% the peak emission in the Ha line has been used 
to diagnose accretion, with values greater than 180-200 
km s~^ indicating accretion for very low-mass stars (e.g. 
White fc Basril 120031: [Muzerolle et all 120031: iNatta et al.l 
2004 iMohantv et all I2005D . Accretion onto low-mass 
stars has also been identified by emission in permitted 
optical lines. Optical forbidden lines are produced by 
outfiows and can therefore be used as indirect accre- 
tion diagnostics be cause outfiows requir e ongoing or re- 
cent accretion (e.g. iHamann et al.l[l99l : iHartigan et al.l 
fl995l) . The presence of excess Balmer and Paschen con- 
tinuum emission is also associated with accretion (e.g. 
IValenti et aDll993t iGullbring eraHll998l ). In this sub- 
section we describe the Balmer jump and excess Balmer 
continuum emission as an accretion diagnostic for the 



individual stars in our sample. In general, the classi- 
fication of stars as accretors or non-accretors based on 
the presence or absence of Balmer continuum emission 
is consistent with previous classifications using alternate 
accretion diagnostics. 

Two mid-M dwarfs in our sample, CFHT 7 and MHO 
7, were previously cla ssified as non-accretors. For CFHT 
7, iGuieu et al.l (|2007l ) measured only weak Ha emission 
and found an absence of excess IR emission shortward 
of and i ncluding the Spi tzer 24/im MIPS bandpass. For 
MHO 7. lMuzerolle et ahi (2003^ measured weak Ha emis- 
sion with a 10% width of 115 km s~^ and no excess 
K-band emission. These two stars have Balmer jumps 
of 0.36 and 0.34, respectively. Both stars are chromo- 
spherically active, with equivalent widths in Ca II H & 
K and He I lines that are larger than is typically mea- 
sured from chr omospheres of older, magnetically-active 
M-dwarfs (e.g. iGiam papa e t all 11981 iGizis et al.l l2002l: 
iRauscher fc Marcvll2 006: AllrMeialj2006D . Theearher- 
type non-accretors LkCa 7 and V819 Tau are bluer and 
have Balmer jumps of 0.44 and 0.46, respectively. 

Many stars in our LRIS sample have much larger ob- 
served Balmer jumps than those non-accretors. The 
stars CM Tau, CIDA 1, 2M0414-h2811, 2M0441-h2534, 
and S0518-I-2327 have Balmer jumps ranging from 2.0- 
3.4 and excess Paschen continuum emission detected at 
> 6000 A. These stars are also associated with strong 
emission in the Balmer series. He I lines, the Ca II IR 
triplet, O I A7773, 8446, and many forbidden lines 

The stars MHO 6, 2M0439+2544, 2M0444-f 2512, 
S0439+2336, 2M0436+2351, and 2M1207-3932 have ex- 
cess Paschen continuum emission detected at 4227 A but 
not at > 5500 A, and five have Balmer jumps between 
0.7-1.4. The two observations of 2M1207-3932 yield out- 
lier Balmer jumps of 3.9 and 5.7. These six targets are 
associated with some strong permitted line emission, in- 
cluding bright Balmer continuum emission and He I line 
equivalent widths of > 1 A, and some forbidden line emis- 
sion. Any Ca II IR triplet emission is faint. 

Each of these stars can be identified as accretors based 
on their high Balmer jump. Similarly, the higher-mass 
stars TW Hya, DF Tau, RU Lupi, and DG Tau have 
Balmer jumps that are well above that of non-accretors 
and can be classified as accretors based on either the 
Balmer jump or the presence of bright emission lines. As 
with other diagnostics, small Balmer jumps with values 
near those of non-accreting templates may be inconclu- 
sive for classifying stars as accretors. V836 Tau has an 
observed Balmer jump of 0.60, which is higher than non- 
accretors with similar spectral type. Many young late 
K and early M stars that are not identified with accre- 
tion show no signs of excess Balmer continuum emission. 
Thus, this excess emission from V836 Tau is identified as 
a clear indication of accretion. 

MHO 5, 2M0455-H3028, CIDA 14, KPNO 11, and V927 
Tau all have observed Balmer jumps consistent with or 
slightly above the non-accretors. Figure 6 shows that 
some excess Balmer emission relative to photospheric 
templates is present from each of these stars, with veil- 
ings at 3600 A that range from 0.14-0.54 (see Table 11). 
Emission line diagnostics suggest that MHO 5, CIDA 14, 
and V927 Tau are actively accreting (see §9.4), while the 
status of 2M0455-)-3028 and KPNO 11 is uncertain from 
other diagnostics. 
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The weak excess B aimer continuum emission from 
these five stars could in principle be produced by the 
stellar chromosphere. Younger M-dwarfs are bluer 
and more chromospherically active than older M-dwarfs 
(jCardini fc Cassatellal l2007f ). which is consistent with 
enhanced Balmer and Paschen continuum emission. 
However, the Balmer jump and elevated Balmer con- 
tinuum emission is not d e tected during stell ar flares 
(jHawlev k, PettersenI 119911 : lEason et alTll992D . which 
may suggest that stellar chromospheres do not produce 
any Balmer continuum emission. In this case, any weak 
excess emission could be attributed to accretion. A larger 
sample of U-band spectra from non-accreting young M- 
dwarfs is needed to definitively rule out a chromospheric 
origin. For the remainder of the paper we treat the ex- 
cess emission from these five stars as upper limits. We 
also suggest the criterion that any observed Balmer jump 
of > 0.5 for a mid-M dwarf should be considered an ac- 
cretor. 

4. MASS ACCRETION RATES 
4.1. Description of Accretion Models 

The accretion emission consists of a strong Balmer con- 
tinuum, a weaker Paschen continuum, and many lines. 
The flux in the accretion continuum is measured for 
the entire 3000-9000 A wavelength range from only the 
strongest three accretors (see Table 4). For the weakest 
accretors the veiling is measurable only shortward of the 
Balmer jump. Calculating the total luminosity of the ac- 
cretion shock requires modelling the observed accretion 
SED to obtain a bolometric correction that accounts for 
unseen emission. 

The H Balmer and Paschen continua and series can can 
be approximated by an isothe rmal, plane-parallel, pure 
hydrogen slab m odel following I Valenti et al.l ()1993[ ) and 
iGullbring et al.l (|T998). These simplistic models were 
developed based on the outdated boundary layer accre- 
tion paradigm but are robust to different geometries for 
the emittin g gas. We use these isothermal plane-parallel 
model from I Valenti et all (|1993f ) because it does not rely 
on a specific morphology of the accretion flow and it pro- 
vides measurements of the accretion lum inosity calcu- 
lated with a methodology c onsist ent with IValenti et al.l 
(jl993) and lCxullbring et all (fl9M). 

The accretion slab models of IValenti et al.l ()1993 f) ap- 
proximate the broadband hydrogen accretion continuum 
with the following free parameters: temperature T, den- 
sity Tie, path length L of the emitting gas (an opacity 
parameter) , turbulent velocity, and filling factor 5 of the 
accretion slab. At low density {n^ < cm~^) the 

Balmer jump is large and determined by the tempera- 
ture of the emitting gas. At higher density stronger H~ 
continuum emission reduces the size of the Balmer jump 
and makes the continuum shape at < 3700 A bluer. The 
Balmer and Paschen continua become redder and the size 
of the Balmer jump increases with decreasing tempera- 
ture, while both continua become bluer and the size of 
the Balmer jump decreases with increasing optical depth 
{Uf. and L). In our fits we vary the temperature, den- 
sity, and path length while setting the turbulent velocity 
to ^ 12 km s~^, whi c h is sm aller than the 150 km s~^ 
used by IValenti et al.l ()1993l ). Larger turbulent veloci- 
ties yield Balmer line profiles that are broader than is 
observed. Substantially smaller turbulent velocities are 



unable to reproduce the shoulder of the Balmer jump, 
located between the real Balmer jump at 3646 A and the 
observed Balmer jump at ~ 3700 A, that is attributed to 
blending of high Balmer lines. The synthetic spectrum 
is then scaled to the measured flux accretion continuum. 
The filling factor of the emission on the stellar surface 
is calculated from the accretion luminosity and temper- 

a ture. 

ICalvet fc Gullbringl (|1998f l developed more advanced 
models of the pre- and post-shock gas at the footpoint 
of the magnetospheric accretion column on the star. In 
their models the Balmer continuum is produced in the 
optically-thin pre-shock gas while the Paschen contin- 
uum is produced in the optically-thick post-shock gas. 
The optically-thin pre-shock gas has an intrinsically large 
Balmer jump. The optically-thick post-shock region has 
a high density and small Balmer jump. The observed 
Balmer jump is therefore not density-sensitive but in- 
stead depends on the relative amount and temperature 
of the pre- and post-shock gas. Our primary goal for 
modelling the measured accretion continuum is estimat- 
ing bolometric corrections to convert the continuum flux 
to an accretion luminosity. ICalvet fc Gullbringl (jl998l ) 
find that the two model variants lead to 2-5% differences 
in the slope of the Balmer continuum and 5-10% differ- 
ences in the slope of the Paschen continua, so that the 
continuum emission from shock mo dels is redder. How- 
ever, our models and the models of ICalvet fc Gullbrind 
()1998f l both successfully explain continuum emission at 
< 5600 A and both underestimate continuum emission 
at > 6000 A (see §4.2). We therefore suggest that any 
differences between the two models in the bolometric cor- 
rection are modest. 

4.2. Application of Models to the Observed Spectra 

We fit the synthetic accretion spectrum plus a pho- 
tospheric template to the unreddened blue spectra be- 
tween 3200-5600 A. The photospheric template is scaled 
to match our estimate for photospheric emission between 
3200-5600 A, which is informed by veiling estimates (see 
Tables 4-5). 

The stars with measurable accretion continuum emis- 
sion longward of the Balmer jump allow us to constrain 
some parameters of the plane parallel slab. Figure 7 
shows that the slope of the Balmer continuum constrains 
T while the Balmer jump constrains n^. Values of T and 
rie for each star are adopted from the center of the accept- 
able contours, with L is subsequently tuned to match the 
strength of the high Balmer lines. Since L and rif. are not 
independent, a lower L may be consistent with a larger 
rif,. In several cases L could not be tuned to yield the 
obs erved flux in the high Balmer lines. In the context of 
the ICalvet fc Gullbrind (|1998f ) shock models, the larger 
Balmer jump may be caused instead by fainter Paschen 
continuum emission from the optically-thick post-shock 
gas relative to the Balmer continuum from the optically- 
thin pre-shock gas. 

The lower Balmer lines are easily affected by emission 
or absorption in the wind and are consistently underesti- 
mated by our models. In magnetospheric accretion mod- 
els of Muzerolle et al. (2001) and Kurosawa et al. (2006), 
emission in the Balmer lines is dominated by the funnel 
flow. If the high Balmer lines are also produed in the ac- 
cretion flow, then the slab optical depths in our models 
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should be larger. The slab optical depth may be smaller 
if instead the wind is optically thick and absorbs a large 
percentage of emission in all of the Balmer series lines. 

The stars S0439+2336, 2M0455+3028, MHO 5, CIDA 
14, KPNO 11, and V836 Tau have weak Balmcr contin- 
uum emission and undetected Paschen continuum emis- 
sion. Since the size of the intrinsic Balmer jump for these 
targets is uncertain, we assume values for rig, T, and I 
based on our other data. We adjust the parameters to 
match the synthetic Balmer series and scale the emission 
so that the spectrum of a photospheric template plus 
accretion matches the total emission spectrum. 

Table 11 lists approximate accretion slab parameters 
for each star, the estimated Balmer jump of the accre- 
tion continuum, the slope of the Balmer continuum, the 
veiling at 3600 and 4000 A measured from these fits, 
and the total accretion slab luminosity. The continuum 
emission from each target is well characterized by tem- 
peratures of - 7000 - 9500 K, densities of (0.3 - 2) x lO^^^ 
cm~^ and lengths of ^ (0.5 — 50) x 10^ cm. In Figure 8 
we show our fits to each spectrum. The synthetic hydro- 
gen continuum is consistent with veiling measurements 
at < 5600 A but underestimates the veiling at > 6000 A 
from 2M0414-t-2811, 2M0441-h2534, CIDA 1, GM Tau, 
and 80518-1-2327 by a factor or 1.2-2. The intrinsic 
Balmer jump of 2M1207-3932 (15 ± 4 and 7 ± 2 for our 
two observations) is particularly large, w hich indicates 
a low H~ opacity. A few o utlier s in the IValenti et al.1 
(|1993f l and iGullbring erall (|1998l ) data also have large 
Balmer jumps. 

4.3. Calculating Mass Accretion Rates and Errors 

The accretion luminosity. Lace, can be converted to the 
mass accretion rate, M, by assuming that the accretion 
energy is reprocessed entirely into the accretion contin- 
uum. The mass accretion rate is then 



M 



1 



Rin 



1.25- 



(1) 



where the factor (1 



1.25 is estimated by as- 



suming the accreting gas falls onto the st ar from the trun- 
catio n radius of the disk, R 

in 

5 i?* ()Gullbring et al.l 
[1991 . The stellar mass and radius are listed in Table 
3 (see also §9.3). We calculate mass accretion rates be- 
tween 2 X 10~^^ to 10~® Mq yr~^ for the very low-mass 
stars and brown dwarfs in our Keck/UilS sample and 
2 X 10-1° - 5 X 10-8 Mq yr-i in the HST/STIS sample 
of higher-mass stars (Table 11). 

The parameters in Equation 1 suffer from systematic 
and random uncertainties related to the data itself, inter- 
pretation of the data, and geometrical assumptions. In 
the following subsections and in Table 12, we describe for 
purposes of explicit clarity these sources of uncertainty 
and how they affect our measurements of A4. 

4.3.1. Errors in Assumed Stellar Properties 

Distance enters into calculations of M (Eqn. 1) be- 
cause Lace oc (P and R (x d, yielding M oc d^. Based on 
the ki nematic distances calculated bv lBertout fc Genoval 
([2006), the distance to any individual Taurus member 
can differ by ~ 15 p c from the standard 140 pc distance 
([Kenvon et al.lll994[ )^. The ~ 11% distance uncertainty 



in Taurus leads to a 0.13 dex in AI. The actual deviation 
in distance may not be random for our Taurus sample be- 
cause the selection of targets here and in previous work 
is influenced by the optical brightness. Targets located 
in front of the cloud will be brighter because of a smaller 
distance modulus and lower extinctions. 

Extinction also affects measures of Lace and R. The 
adopted 0.5 mag uncertainty in Ay (see §9.2) leads to a 
0.32 dex uncertainty in Lace and a 0.08 dex uncertainty 
in R for a total of 0.4 dex uncertainty in M. If Ry > 3.1, 
then the extinction curve would be flatter and our Ay 
would be underestimated. 

The M measurements also depend on the ratio 
which is calculated from the effective temperature and 
theoretical evolutionary tracks (§9.3). The tempera- 
ture is uncertain by '--^ 75 K, including both the un- 
certainty of roughly 0.5 spectral type subclasses and 
the conversion of spectral type to temperature (§9.1). 
In the Baraffe et al. (1998) pre-main sequence evolu- 
tionary tracks at 2 Myr and with temperatures from 
iLuhman et al.l ([2003f ). a difference of 0.5 spectral classes 
corresponds to 25% larger mass and a 10% smaller radius. 
Alternate pre-main sequence evolutionary tracks from 
iD'Antona fc Mazzitellj (|1998D would increase our mass 
estimates 1.1 -1.7 times larger than the [Baraffe et al.l 
(1998) tracks. I White fc Hillenbrand! (|2004[ ) suggest that 
evolutionary tracks underestimate stellar masses by 30- 
50%. 

All sources except V927 Tau and DF Tau are assumed 
to be single. If any sources are binaries, the lower lumi- 
nosity per star would imply that we have overestimated 
R. The shallower potential well would result in AI being 
1.4 times lower than is calculated here. In cases where 
both sources are accretors, the accretion luminosity and 
mass accretion rate is the total rate onto both stars. 

4.3.2. Uncertainty in Bolometric Corrections 

Converting the measured accretion flux into an ac- 
cretion luminosity requires a bolometric correction for 
unseen emission (§4.2). However, the shape of the 
broadband accretion continuum is not well under- 
stood. Our single-temperature slab models underes- 
timate the veiling at 6000 — 9000 A by factors of 
1.2-2, as whic h also occurs for the multi-temperature 
shock models ([Calvet fc Gullbrin e 1998) . The veilin; 
at even longer wavelengt hs (e.g.. White fc Hillenbrani 
12004 [Edwards et al.|[200l) also seems larger than can be 
explained by the isothermal slab that fits well the Balmer 
continuum and short-wavelength Paschen continuum. 

If we hypothetically double the accretion continuum 
flux at A > 6000 A, the total accretion luminosity in- 
creases by 0.15 dex. If instead the flux per A is constant 
between 6000-20000 A (despite the presumed presence 
of the Paschen jump), then the accretion luminosity in- 
creases by 0.3 dex. A smaller uncertainty of ^ 0.1 dex 
in Lace is introduced by the uncertainty in temperature 
(~ 1000 - 2000 K) of the plane-parallel slab. 

4.3.3. Exclusion of line emission 

Our M estimates rely on accurate calculations of the 
total accretion energy that escapes from the shock. Such 
emission is seen in both continuum and line emission, but 
we include only the total continuum emission from the 
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accretion slab in estimating Lace- The observed Balmer 
line fluxes (Table 7) total 0.2 times Lace for most of the 
accretors but are roughly equivalent to Lace for the M7- 
M8 accretors (see Tables 7 and 11). The formation of 
and flux in lower Balmer lines is complicated by emis- 
sion and absorption in the accretion funnel flow and in 
outflows. In principle other lines (e.g., Ca II IR triplet, 
Lyman series, FUV and X-ray lines) may also need to be 
included. For TW Hya, X-ray and UV lines account for 
30 — 45% of the total accretion flux, most of which is 
in H I Lya, but most have fluxes t hat when summed are 
likely < 20% of the accretion flux (jKastner et al. 1 120021 : 
iHerczeg et al.|[200^ . Our exclusion of line emission from 
the accretion luminosity is consistent with previous Lace 
estimates and is essential for simplifying these calcula- 
tions and comparisons with existing M estimates. 

4.3.4. Geometric assumptions 

In the magnetospheric accretion paradigm, gas ac- 
cretes from the inner disk truncation radius. Rim onto 
the star along magnetic fleld lines, where Rin is deter- 
mined by where the m agnetic field intercepts the disk. 
iGullbring et al.l (|1998[ ) noted that infall energy depends 
on Rin and therefore included the factor of (1 — -^) in 
Equation 1. The Rin ranges from 3-10 i?* b ased on mag- 
net ic field strengths of higher-mass CTTSs (iJohns-KrullI 
[200 Tj) . We assume Rin = ^R* for all targets in this sam- 
ple for consistency with the same assumption made by 
iGuUbring et al.l (| 19981 ). A factor of 2 error in leads 
to a ~ 0.05 — 0.12 dex error in the mass accretion rate. 

As the gas accretes along the magnetic field lines 
and shocks at or near the stellar photosphere, some of 
the emission fr om the hot sp ot will be dir ected at the 
star. W e follow I Valenti et all (|1993l) and IGullbring et all 
(|1998f l in assuming that all of the accretion energy es- 
capes from the star as hydrogen continuum emission. 
By this assumption, any accretion emission directed at 
the star is either coherently re-radiated by the star or 
provides some heating to the nearby g as, which then 
cools by h ydrogen continuum emissio n. iHartigan et al.l 
(|1995f ) and IHartigan fc Kenvonl ()2003l ) instead assumed 
that half of the Balmer and Paschen continuum emission 
is absorbed by the star, which increases the calculated 
M by a factor of two. The accretion luminosity may also 
not be isotropic. 

The accretion hot spot is assumed to be the source 
for all Balmer and Paschen continuum emission from 
our sample. No excess emission is seen in a sample of 
non-accreting K-dwarfs that serve as templates for the 
higher-mass accretors discussed here. Thus, the weak ex- 
cess Balmer emission from V836 Tau that is attributable 
to accretion can be accurately measured. Whether the 
faint excess Balmer emission from 2M0455-I-3028, CIDA 
14, MHO 5, KPNO 11, and V927 Tau is produced by 
accretion or may be chromospheric is uncertain. Using 
V927 Tau as a template may lead us to underestimate the 
upper hmit of M onto 2M0455-H3028 by 0.2 dex. Weak 
Balmer continuum emission in the template spectrum 
has a negligible effect for targets with larger Lace/Lboi- 

4.3.5. Cumulative Effect of Uncertainties on Estimates of 
Mass Accretion Rates 



The uncertainty in bolometric correction, exclusion of 
emission lines, and the re-radiation of H Paschen and 
Balmer continuum emission by the star can only increase 
the estimated accretion rate, while the uncertainty in bi- 
narity and pre-main sequence tracks can only decrease 
the estimated accretion rate. When combined, these er- 
rors suggest that our accretion rates have a relative un- 
certainty of ~ 0.6 dex and underestimate the accretion 
rate by a factor of 1.2 (or 2.4, if one assumes no re- 
radiation of the accretion luminosity by the central star) . 

This error analysis is not a rigorous assessment of the 
uncertainty in M at a specific confidence interval. The 
uncertainty differs for various targets depending on the 
accuracy of stellar parameters, and for many targets the 
error is dominated by uncertainty in Ay. For 2M1207- 
3932, which has an accurate Ay = 0, the M is likely un- 
derestimated by a factor of 2 (or 4, if the accretion emis- 
sion is not re-radiated), primarily because we exclude line 
emission in calculating M . The standard deviation of 

0.3 dex about this factor of 2 is dominated by system- 
atic uncertainties in the bolometric correction. System- 
atic uncertainites of this magnitude are also present in 
all previous UV-excess or optical veiling estimates of M . 
Several uncertainties, including in the extinction mea- 
surements, could be reduced with the addition of near- 
UV spectra, better photospheric templates across the en- 
tire wavelength region, a better understanding of veiling 
at long wavelengths. 

Intrinsic variability in accre tion rate can also a ffect 
non-simultaneous obserations. I Valenti et al.l ()2003f ) de- 
scribe variability of near-UV emission in repeated lUE 
spectra of young stars. The average near-UV variabil- 
ity is 29% for the 15 K1-M2 accretors that were ob- 
served at least 5 times with lUE. The near-UV emission 
from these stars is dominated b y and should correlat e 
directly with the accretion rate (|Gullbring et al.ll2000f) . 
Variability of 0.1-0.2 dex is also commonly seen in re- 
peated veiling measurements of DF Tau, TW Hya, and 
AA Tau ( Johns-KruU fc BasTill 19971: iBatalha et al.ll2002l : 
iBouvier et al.l I2007D . IScholz fc Javawardhanal (|2006f ) 
confirm that variability in accretion rate is also seen for 
brown dwarfs. We infer that the accretion rate is variable 
by - 0.2 dex. 

5. RELATIONSHIPS BETWEEN LINE AND CONTINUUM 
ACCRETION DIAGNOSTICS 

The Balmer continuum can be difficult to detect be- 
cause of large extinctions or instrument limitations. 
Other accretion diagnostics are therefore usually more 
accessible and can be used to infer accretion rates. In 
the following subsections we derive and evaluate empiri- 
cal relationships between emission lines and accretion lu- 
minosities and rates. In §5.1 we calculate accretion rates 
from published data to incorporate additional objects in 
our analysis in §5.2-5.4. We then compare the accretion 
luminosity and resulting accretion rate with Ha line pro- 
file measures of accretion (§5.2), line luminosities (§5.3), 
and line surface fluxes (§5.4). 

5.1. Calculating accretion luminosity from optical 
veiling measurements 

Although ours are the flrst Balmer jump measurements 
of young brown dwarfs, detections of excess Paschen 
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continuum emission have been published. Calculating 
M from these objects increase the overlap of M mea- 
sured using excess co ntinuum em ission and Ha line pro- 
file modelling (§5.2). E raus et a l. (2006) use broadband 
HST photometry to measure excess V-band emission 
from several young brown dwarfs. The V-band excess is 
significant when ry > 1-3, which restricts usage to high 
mass accretion rates relative to the given stellar mass. 
We use log iacc = — 29.59-f logF545o (calculated from Ta- 
bles 4 and 11) to measure M of 10^^ — 10^""^^ Mq yr~^ for 
GM Tau, KPNO-Tau 4, KPNO-Tau 6, and KPNO-Tau 
12 (Table 13). The M from GM Tau based on V-band 
excess is nearly identical to that measured more directly 
from our low-resolution optical spectra. However, the 
true M was likely 20% larger during the HST observa- 
tion because the V-band excess was ~ 0.2 mag brighter 
t han was estimated here. 

iHartiean fc KenvonI (l2003h analyzed HST/STIS spec- 
tra of binaries with early M spectral types and masses 
mostly within the ^ 0.25 — 0.55 Mq range, includ- 
ing veiling measurements at 6100 A and luminosities 
in [O I] A6300, Ha, and Ca II A8542 fines. Using 
log Lace = — 29.67-|-logF6i9o calculated empirically from 
our Keck //LRIS spectra (see Tables 4 and 11), we 
find Lnr.r. that are . 28 — .75 times those calculated by 
iHartigan &: KenvonI (|2003f ). We adopt these lower Lace 
for further analysis. 

5.2. Comparing UV-excess Accretion Rates to Ha Line 
Profile Models 

Strong Ha emission is usually attributed to accre- 
tion. At high-resolution the line profile often in- 
cludes P-Cygni, inverse P-Cygni, or self-reversed ab- 
sorp tion that corrupts equivalent width measu rements 
(e.g. iReipurth et al.|[l99l iMuzerolle et al.ll2003D . Muze- 
rolle et al. (1998, 2001) appfied magnetospheric accretion 
models to calculate the formation and radiative transfer 
of Ha lines. These models treat the line profile, scaled 
by stellar parameters, as a proxy for M. MuzeroUe et 
al. (2003, 2005) applied these models to very low-mass 
stars and brown dwarfs to estimate accretion rates. Fig- 
ure 9a and Table 14 compare mass accretion rates for 
stars with M measured from both Ha line profile models 
and UV-excess measurements. 

The UV excess accretion rate measurements are ~ 0.84 
dex larger than the Ha accretion rate measurements with 
a standard deviation of 0.3 dex, if MHO 5 and CIDA 14 
are ignored. If both the Ha and weak excess Balmer 
emission from those two stars are attributable to accre- 
tion, then the UV excess measures are 0.6 dex larger 
than the Ha accretion measures but with a larger (0.5 
dex) source-to-source variation. 

This difference is consistent with the systematic 
methodological differences between the two methods. 
However, that the systematic errors are in opposite direc- 
tions for the two methods is unlikely. For the UV excess 
method several uncertainties (see §4.3) would serve to 
increase our M estimates and therefore the discrepancy 
between the two methods. The uncertainty in the Ha M 
calculations of ^ 3 — 5 are dominated by uncertainty in 
the disk trun cation radius of the in ner disk, assumed to 
be 2.2-3 ii* (|Muzerolle et al.ll2003[ ). The model Ha line 



width is proportional to M^^'^, so overestimating M^ will 
lead to underestimating M for both methods. Intrinsic 
variability in M can explain much of the scatter once 
the systematic variation between these two methods is 
corrected. 

In the MuzeroUe et al. Ha line profile models, the ac- 
cretion rate is estimated from the line profile along with 
the stellar mass and disk inclination. Although M is cor- 
related with the Ha line width in the models, different 
M may vary by more than an order of magnitude fo r 
equivalent 1 0% wi dths (see also iKurosawa et "al]|2006f) . 
iNa tta et al] (|2004l ) found an empirical relationship be- 
tween M and Ha 10% width, which is used in the ab- 
sence of Ha line profile modelling. We find a 0.86 dex 
scatter between the M mea sured here from UV excess 
and M from the lNatta et all (2004) Ha relationship (Fig. 
9b) and a similar amount of scatter for a similar best-fit 
fine'' from Table 15. 2M1207-3932 is a particular out- 
lier on this plot, with an M estimated from the Ha 10% 
width that is two orders of magnitude larger than that 
measured from the UV excess. 

5.3. Empirical relationships between line and continuum 
luminosities 

In §3.4 we qualitatively describe that excess continuum 
and line emission are produced by accretion processes. 
Lines from accretors exhibit a wide range of profiles 
in high-resolution spectra, with eq uivalent widths that 
are co rrelated with v eiling ( e.g. iHamann fc PerssonI 
1992t iBat alha et al.' '199 6t IMuzerolle et al.l 119981: 
Beristain e t al. 2001; Moh antv et al.l I2005D . Figure 
10 shows that the hmrinosity in many lines {Lune) is 
tightly correlated with accretion continuum luminosity, 
although this relationship does not require that the 
emission lines and continuum are produced by the same 
gas. We improve the fits for the Ha, [O I] A6300, and 
Ca II A854 2 lines by suppl ementing our data with results 
from Hartigan fc KenvonI (j2003i) . All line and accretion 
luminosities used in this subsection were measured from 
simultaneous data. 

Linear fits of the relationship between log Lace and log 
Liine are calculated using the statistical package asurv 
(jFeigelson fc NelsonI IT985I) . with upper limits ignored 
(see Table 15 and Figure 10). Most upper limits are 
not significant and their inclusion does not substantially 
change the results. 2M0455-h3028, KPNO 11, CIDA 14, 
MHO 5, or V927 Tau (squares in Figure 10) are not 
included in these fits because the source of the excess 
Balmer emission is uncertain. For several lines at > 5600 
A, the high M range for these fits is dominated by CIDA 
1 and GM Tau, both of which have large uncertainty in 
Ay and therefore M. The fits based on few data points 
should be used with caution. 

The Ha-Lace correlation has more scatter than the cor- 
relation between higher Balmer lines and Lace because 
the larger oscillator strength of Ha can lead to stronger 
wind absorption. The Balmer and Ca II H fc K line lu- 

* We find that M = -14.0-|-0.014(Ho 10% width) from the data 
in Table 14, excluding upper limits, and M = -12.6 -I- 0.010(Ha 
10% width) if we combine this data with additional M calculated 
from the Ca II A8662 line fluxes from Mohanty et al. (2005, see 
also §5.4). 
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minosities from KPNO 11, 2M0455-I-3028, MHO 5, and 
CIDA 14 lie above the linear fit to accretion luminosity. 
Since these four stars are tentative weak accretors, such 
line emission is likely produced by the chromosphere. 
The relationship between the Balmcr and Ca II K lines 
and the accretion luminosity d erived from our data i s 
consistent with the results from iGullbring et al.l (|1998D . 
The luminositie s in these lines are systematically ~0.5 
dex lower in the IValenti et al.l (|199S) sample than in the 
IGullbring et al.l ( 19981 ) sample.^ The nature of this dis- 
crepancy is not known. 

A large amount of scatter in the relationship between 
[O I] line luminosity and accretion luminosity arises be- 
cause the [O I] e mission line forms in an accretion- 
jowered outflow (iH amann et al.l |1994; Hartiga n et al 
The comparison of the iHarti gan fc Kenvon 



sample and our sample may be corrupted be- 
cause [O I] emission i s typically spatially-extended (e.g. 
IHartigan et al. I l2004f ). Any extended [O I] emission 
could lie outside the narrow 0'.'2 slit width used by 
(|Hartigan et al.llT995f ). but our I'.'O slit width includes 
any forbidden line emission related to the star. The 
outflow line luminosities are correlated with accretion 
rate, though strong accretors sometimes drive weak out- 
flows while some weak accretors may drive powerful out- 
flows. For example, the outflow lines from GM Tau and 
2M0414+2811 are substantially stronger than those from 
CIDA 1, despite the higher mass accretion rate onto 
CIDA 1. Similarly, many bright outflow lines are de- 
tected from 2M0444-t-2512 despite a modest mass ac- 
cretion rate. The He I and Na I D li nes are usually 
associated with both infal l and outflow (IHartigan et al.l 
ll995HBeristain et al.ll200lHEdwards et al.ll2006l ) but here 
are more tightly correlated with the accretion luminosity 
than [O I]. 



5.4. Comparison with M measurements from emission 
line fluxes 

In §5.3 we derive empirical relationships between 

line and accretion conti nuum luminosit y. Lacking 

flux-calibrated spectr a, iMuzeroUe et al.l (|1998f ) and 
iMohantv et al.l ()2005( ) converted the Ca II equivalent 
width into line flux per stellar surface area and then cor- 
related the line flux directly with M. Such a correlation 
naturally corrects for radius uncertainties introduced by 
errors in extinction or distance. Compared with equa- 
tion 1, correlating ^'^S^ with M will be in error in M by 

a factor of 

Figure 11 compares the Ca II A8542, Ca II A8662, and 
He I A5876 line fluxes with accretion rate for the low- 
mass accretors presented here, high-mass accretors from 
MuzeroUe et al. (1998 for the A8542), and low- and high- 
mass accretors from Mohanty et al. (2005 for the A8662 
line)^. Those data have M and Ca II line fluxes measured 
from non-contemporaneous observations and will suffer 
from variability. The i^(Ca II A8542) - M and i^(He I 



H 



^ In IValenti et al.l I I1993I) and here, the Balmer and Ca II 
K lines are measured directly from the spectrum while 



A5876) — M relationships are shallower than was previ- 
ously measured from only the high-mass sample, while 
the J^(Ca II A8662) — M relationship is similar to that 
derived by Mohanty et al.(2005). The standard devia- 
tions between M predicted from these line correlations 
and the M measured from UV continuum excess are 1.0, 
0.71, and 1.1 dex for Ca II A8542, Ca II A8662, and He I 
A5876 Hue fluxes, respectively. 

6. DISCUSSION 

6.1. Selecting methods to measure M 

In the magnetospheric accretion paradigm, the 
Paschen continuum emission is produced in optically- 
thick post-shock gas, the Balmer continuum emission 
is produced in optically-thin pre-shock gas, and the 
Balmer line emission is mostly produ ced primarily in 
the magnetosph eric accret ion co l umn (iHartmann et all 
1991 ICalvetfc GuUbrii^i [l998l : IMuzerolle et al.l 120011 : 



Kurosawa et al.li200 6). In the simplistic plane-parallel 



IGullbring etaLl i ITQQI) measures line fluxes after the subtraction 
of a photospheric template. These two different methods would 
slightly exacerbate the discrepancy between the line luminosities. 
6 We calculate logAf = 1.28 log F{Ca II A8542) - 16.6, logM = 



models all three accretion diagnostics are produced by 
the same gas. In either case, the accretion energy is re- 
processed into Balmer and Paschen continuum emission, 
which can be measured and subsequently converted into 
accretion rate. Previously, the accretion continuum mea- 
suren ients were a pplie d only to accretors with ^ 0.3 — 2 
Mq. iMuzerollTet al.l (|2003l ) developed models to mea- 
sure M from Ha line profiles for only the very low-mass 
stars and brown dwarfs. At higher mass accretion rates 
the Ha emission line profile is difficult to use for estimat- 
ing M because emission and abso rption in strong stellar 
outfiows corrupt the line profile (jMuzerolle et al.l [20011 : 
Kurosawa et al.l [20061 ). We provide the first measure- 
ments of the Balmer continuum emission from very low- 
mass stars and brown dwarfs and describe this emission 
as qualitatively similar to that from higher-mass accre- 
tors. The UV excess measures of M are ^4 — 7 times 
larger than those from Ha line profile modelling with 
some source-to-sourcc differences that may be method- 
ological or related to intrinsic source variability. 

In the context of magnetospheric accretion models, 
that the excess accretion continuum and the Ha emis- 
sion yield different M may not be surprising because they 
form in different locations. The primary advantage of the 
Ha models is that the line profile is independent of the 
intrinsic line luminosity and onl y moderately dependen t 
on the stellar mass. According to IMuzerolle et al.l (|2003f) . 
the uncertainty in M from Ha modelling is a factor of 
3-5 and is dominated by uncertainty in the size of the 
magnetosphere. The factor of 4 random error in UV 
excess measures of accretion is dominated by errors in 
extinction and distance. The UV excess method relies 
predominantly on directly observable parameters, with 
modelling required only for calculating a moderate bolo- 
metric correction to derive Lace- The accretion energy 
is converted to M by equation 1. In contrast, accretion 

1.03 log F(Ca II A8662) - 15.2, and logAf = 1.63 log F(Hc I 
A5876) — 17.9 from line fluxes and UV excess measures of accretion. 
We also note that the relative fluxes in the three Ca II IR triplet 
lines are usually consistent with the optically-thick limit of 1:1:1 
l|Herbig fc SoderblomllQSa iHamann fc Persson'1989'). whether the 
emission is in a broad or narrow line (Batalha ct al. 1996). Even at 
high-resolution all three Ca II IR triplet lines blend with Paschen 
lines. 
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rates from the Ha line profile rely on a correct interpre- 
tation and modelling of the line profile in the magneto- 
spheric funnel flow, and may be complicated by outflows, 
emission from hotspots, complex magnetic fields, or un- 
expected geometries. 

Figures 12-13 compares the spectral properties of the 
excess hydrogen continuum emission versus stellar mass 
and mass accretion rate. The size of the Balmer jump 
tends to increase with decreasing M (and M), consistent 
with lower ne (~ 3 x 10^3 cm'^ for 2M1207-3932 com- 
pared with ~ 2 X 10^'* cm~3 for higher mass accretors) 
in single-temperature slab models or with a larger con- 
tribution of emission from the optically-thin pre-shock 
gas relative to optically-thick post-shock gas in the mag- 
netospheric accretion models. The slope of the Balmer 
continuum, which may be affected by our atmospheric 
correction, tends to increase with decreasing M (and M). 
This is consistent with a lower temperature (~ 7500 K 
for 2M1207-3932, compared with - 8000 - 12000 K for 
higher- mass accretors) in the plane-parallel slab models. 
MuzeroUe et al. (2003) describes that the temperature 
in the magnetospheric accretion flow must be at least 
10^ K to have sufficient opacity in Ha to produce the 
broad line profiles. The percent of the stellar surface 
subtended by the accretion hot spot (5, see Table 11)^ is 
typically 0.01-1% for low-mass stars and brown dwarfs, 
relati ve to 0.1-10% for higher- mass stars (jValenti et al.l 
Il993l ). These results should be fully accounted for as 
tests of the magnetospheric accretion models. 

For regions with low extinction the UV excess method 
is the most accurate measure of accretion luminosity. 
The accretion luminosity can also be measured from the 
Paschen continuum at 3700-8000 A, albeit with larger 
uncertainties and reduced sensitivity to small accretion 
rates. UV excess measures of accretion can be applied to 
even lower mass objects than studied here, limited only 
by telescopic sensitivity. Ha line modelling may be a 
sufficient substitute for low-mass stars and brown dwarfs 
when UV excess measures are not possible, either due to 
large or uncertain extinction or distance, or instrumen- 
tation limitations. The Ha models at present cannot be 
applied to large M due to complex line profiles. They 
also cannot be applied to M < 10~^^ Mq yr~^, below 
which the opacity in the Ha line wings may be insuffi- 
cient to measure accretion (|Muzerolle et al.ll2003l) . How- 
ever, Ha models may be better suited than UV excess 
measurements for small M onto higher mass/luminosity 
accretors, where the contrast between the Balmer con- 
tinuum and photospheric emission is small. For exam- 
ple, the excess Balmer emission is difficult to detect onto 
the objects CIDA-14, V927 Tau, and 2M0455-h3028, but 
the 10% width of the Ha line is large enough to suggest 
accretion. Additional comparisons of M between these 
two methods are required to better characterize the dif- 
ferences. 

Secondary diagnostics, such as the Ca II IR triplet, 
Pa/3, or Br7 lines, are related to Lace or, less accurately, 
M by empirical relationships and can provide useful 

The estimates of S are inversely correlated with L, the length 
of the accretion slab, and T, the temperature of the accretion slab. 
The calculated L may be non-p hysical, and in this w ork are larger 
here than the ~ 10^ typical for lValenti et afl l)1993I V As a result, 
the difference in 5 between the two works may be underestimated. 



M es t imates (see also iMuzeroUe et"aLlll998l : iNatta et al.l 
12004 iMohantv et~al] l2005l ). The manv random uncer- 
tainties inherent in both the UV excess and Ha methods 
are reduced if the sample size is large enough. However, 
systematic uncertainties in the both accretion measures 
are propogated into these correlations. Some uncertainty 
is also introduced into these relationships because the 
line emission may be produced in different gas than the 
accretion continuum and because wind absorption can 
reduce the line flux. Of the line luminosity-accretion lu- 
minosity relationships analyzed in Table 15, the accre- 
tion luminosities estimated from the high Balmer lines 
and the He I A5016 line have the least amount scat- 
ter. The line luminosity-accretion luminosity relation- 
ships are preferable to line flux-Af relationships (Fig. 
11), which have scatters of ^ 0.75 — 1.0 dex. That this 
amount of scatter is similar to the scatter in the M — M 
relationship (see §6.2) suggests that the line flux-M re- 
lationships are not accurate. 

Many of the commonly used lines (Ca II IR triplet, 
Pa/3, Br7) are less sensitive to M/M than either the 
Ha line profile or excess UV emission. They therefore 
provide a less complete picture of the full range of ac- 
cretion. The completeness depends on mass if detection 
limit of each diagnostic varies with spectral type. These 
diagnostics are often more accessible than modelling Ha 
profiles or measuring the UV excess and are particularly 
useful for va riabilit y studies of a individ ual objects (e.g. 
lAlencar eFal.„2005i:IAzevedo et al.ll2006[) or in large sam- 
ples (|Natta et al.ll2006f) " 



6.2. Accretion rate versus mass 

Figure 14 shows M versus M for Tauru s accre - 
tors with Ay < 2. 5 fro m iGullbring erahl ^30), 
Hartigan fc KenvonI (l200l. iCalvet et afl ^04), 
Muzerolle et al. I (|2005D . and from the low-resolution 
Keck I and HST spectra analyze d in this work*. We use 
asurv ([Feigelson fc Nelsonlll985[ ) to calculate that 

logM = (1.87 ± 0.26) logM- (7.73 ±0.17), (2) 

which is similar to previous estimates of M (x AI^-^^'^-^ 
(|Muzerolle et al.ll2005t IMohantv et"al]|2005t [Natta et all 
120061 ). The standard deviation of the fit to M is 0.8 
dex. This standard deviation is too large to be explained 
by intrinsic accretion variability or random methodolog- 
ical uncertainties, which suggests that stellar parameters 
other than mass also affect the accretion rate. A lack of 
completeness in our sampling of the true range in mass 
accretion rate for a given mass can bias the M — M" 
relationship by target selection and by sensitivity limits. 

We specifically selected several targets (CIDA 1, GM 
Tau, and 2M0414-I-2811) based on our expectation of 
high M. Most other stars were selected based in part 
on visual magnitude, which may tend to yield sources 
with higher M. If this selection bias is only present in 
the current data, which focuses on the low-mass range, 
then the a in the M — M" relationship will be underesti- 
mated. However, a similar bias could also be present in 



* The IHartigan &: KenvonI l l200l) data are used to re-calculate 
accretion rates, as in §5.1. The IMuzeroUe et al.l I I2005I) accretion 
rates are used only if they are not duplicates from this work and 
are multuplied by a factor of 5 (see §5.3). Duplicate measurements 
of M for the same star are averaged. 
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previous M measurements from samples of higher-mass 
stars. 

The M from higher-mass accretors may be incomplete 
because sensitivity to |y improves to smaller masses for 
both UV excess and Ha accretion measures. Figure 15 
shows that accretion luminosity and the measured 
decreases with mass. The improved sensitivity to low 
at smaller masses results because the percentage of 

total photospheric emission that escapes at < 3700 A is 
m uch less for mid-late M-d warfs than for late-K dw arfs. 
In iValenti eFaLl (|1993D and iGullbring et all (|1998D . the 
minimum value of is 0.008 and 0.018, respectively, 
with most targets well above that limit. This detection 
hmit leads to ^ 5 - 10 x 10" Mq yr^^ for a 2 Myr 
old 0.7 Mq star with L = 0.57 Lq. For the lower-mass 
accretors studied here, accretion from 2M0455-I-3028 is 
at our detection limit with f""^"" '-^ 0.001. For a 2 Myr 
old, 0.1 Mq brown dwarf, the lowest detectable mass 

accretion rate is ~ 10"^^ Mq yr"^. The ^ ratio is 
therefore about 5.5 times lower for a 0.1 Mq accretor 
than for a 0.7 Mq accretor. Incompleteness at higher 
masses would cause us to overestimate a in the M — M" 
relationship. On the other hand, inclusion of the five 
tentative accretors would increase a to 2.08. 

7. CONCLUSIONS 

We analyzed blue spectra of 18 accreting very low-mass 
stars and brown dwarfs obtained with LRIS on Keck I. 
Ours are the first data on the Balmer continuum and 
high Balmer line emission in this mass range. Most of 
this sample was selected based on the presence of ac- 
cretion identified with other diagnostics. Our observa- 
tions are compared with archival HST/STIS spectra of 
the Balmer continuum from several higher-mass stellar 
accretors. We obtain the following results: 

1. We detect excess hydrogen Balmer continuum emis- 
sion from 16 of the 18 mid-late M-dwarfs in our sample. 
11 of these objects also show hydrogen Paschen contin- 
uum emission. This continuum is attributed to accre- 
tion for most of our objects. Weak Balmer continuum 
emission onto earlier-type stars, including the K6 dwarf 
V836 Tau, can be attributed to accretion because non- 
accretors show no evidence for chromospheric Balmer 
continuum emission. However, we cannot rule out a chro- 
mospheric origin for weak excess Balmer emission from 
mid-M dwarfs. 

2. Measuremeiits of the Balmer and Paschen continua 
emission yield M — 10"^^ — 10~^ for very-low mass 
stars and brown dwarfs with masses from 0.024-0.17 Mq. 
Most of these accretion rates are smaller than those typi- 
cally measured for accreting ~ 0.3 — 1.0 Mq stars. Anal- 
ysis of the random and systematic errors indicate that 
these M may be underestimated by a factor of 1.2 (or 
2.4 if the hydrogen continuum emission incident upon 
the star is reprocessed into photospheric emission), with 
a relative uncertainty of ^ 0.6 dex. Comparing these 
M with M from other studies requires accounting for 
different assumptions regarding radiative transfer. 

3. The Balmer jump tends to be larger and the Balmer 
continuum tends to be redder for stars with lower mass 
accretion rates. These properties indicate lower opacity 
(with log rig = 13.5 for 2M1207, compared with logrie ~ 



14.3 for stars with higher M in the plane-parallel slab 
models) and lower temperature (with T ^ 7000 K for 
2M1207, compared with ~ 9500 for stars with higher 
M). 

4. The M calculated here are well correlated with but 
~ 4— 7 times larger than the M calculated from published 
models of the Ha line profiles. We suggest that this offset 
is methodological, although the difference is within the 
range of the combined errors for both methods. The M 
measured here also correlates with Ha 10% width but 
introduces 0.9 dex uncertainty in M predicted from Ha 
10% width alone. 

5. We provide empirical relationships between Lace 
and Liine, which are often well correlated. In particular, 
estimating Lace from the luminosity in the high Balmer 
and many He I lines introduces an uncertainty of only 
0.2-0.35 dex. On the other hand, estimating the Lace 
from line surface ffuxes alone introduces an uncertainty 
of ~ 0.7 — 1.1 dex in M estimates. 

6. We find that logM oc M^-^"^ for Taurus stars with 
Ay < 2.5, which is consistent with previous estimates. 
This relationship may be biased if the higher sensitivity 
of M/M at lower M causes the percentage of accretors 
that are detected to also vary with M. 
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9. APPENDIX 

9.1. Photospheric Emission 

The photospheric emission from M-dwarfs is domi- 
nated by strong TiO and CaH bands at A > 4800 
A that can be used to estimate spectral type (e.g. 
iKirkpatrick et al.lll993f l. Table 3 lists spectral types for 
our ta rgets estimated from the TiO 7140 and 8465 A in- 
dices (jSlesnick et al.|[2"006[) and by visual comparison to 
existing Keck //LRIS spectra of M4.5-M9.5 M-dwarfs^. 
Most spectral types are consistent with existing litera- 
ture estimates. We classify V927 Tau as an M5 star, 
compared with previous classifications of M3, M4.75, 
and M5.5. 2M0455-f 3028, S0439-I-2336, and MHO 6 
have the same spectral type as V927 Tau and are clas- 
sified as M5. MHO 7 is a slightly later spectral type 
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than V927 T au and is clas si fied a s M5.5, 0.25 classes 
later than in iBriceno et alJ ()2002l ). The photospheric 
emission from CIDA 14 is nearly identical to MHO 7. 
S0518+2327 is about the same as MHO 7 and classi- 
fied as M5.5, while KPNO Tau 11 shghtly later than 
MHO 7 and is classified as an M5.75. MHO 5 is between 
MHO 7 and CFHT 7 and is cl a ssified as M6. We confirm 
spectral types from iLuhm^ (|2004D for 2M0438+2611, 
2M0 439+2 5 44, 2M 0441+2534, and 2M0444+2512, and 
from i GizisI ()2002D for 2M1 207-3932. We conservatively 
assign a 0.5 subclass uncertainty to all spectral types. 

The spectral types for GM Tau, CIDA 1, and 
2M04144-2811 are difficult to determine because of high 
veiling at A > 7000 A. In order to make some estimate 
we assume that the accretion continuum emission is con- 
stant in flux per A across short wavelength intervals. The 
depth of absorption bands are then compared to spec- 
tral type standards to estimate the spectral classification. 
The photospheric emission of CIDA 1 is consistent with 
M5, GM Tau with M5.5, and 2M0414-f2811 with CFHT 
7. The spectral types of CIDA 1 and GM Tau are 0.5 
and 1 class ear l ier, re spectively, than that measured by 
iWhite fc BasTil (l2003h and are both assigned an error of 
1.0 subclasses. 

9.2. Extinction 

Extinction estimates are important to convert the 
measured flux in the accretion continuum into a lu- 
minosity. Many of our targets were selected for ob- 
based on low extinction. We adopt ex- 
of ~ 0.4, 0.1, 0.0, and 0.0 to the non- 
V927 Tau, M HO 7, CFHT 7, and the 

2M 0444+2512 jK envon fc Hartmann 1995'; 

IBriceno et all l2002l : iGuieu et a l. 2007; Luhman 2004). 
We confirm these estimates with an accuracy of ^ 1.0 
mag by comparing our flux-calibrated observations to ex- 
isting Keck //LRIS red spectra (6000-10000 A) of pho- 
tospheric templates. 

The photospheric emission from the rest of our tar- 
gets are then compared to those four templates to cal- 
culate an extinction. For accretors, we also add a flat 
accretion continuu m at red wave l engths . We use the ex- 
tinction law from Car delli et al.l ()1989l ) and a total-to- 
selective extinction {Rv) of 3.1. Some evidence suggests 
that Rv > 3.1 in the Taurus Molecular Cloud because 
grains in the molecular cloud may be larger than thos e 
in the average interstellar medium (jWhittet et al.lf2"001f ) . 
If Rv is large, the slope of extinction versus wavelength 
is shallow and Av will be underestimated. The spec- 
trum of 2M0438+2611 is likely seen through an edge- 
on disk and suffers f rom gray extinction and scattering 
(jLuhman et al.|[2007l ) that is not corrected for here. 

We conflrm our extinction estimates for sources with 
strong veiling by assuming that the Paschen continuum 
between 4000-6000 A is flat (see §3.1), which is roughly 
consi stent with observations of the higher- mass accretors 
(e.g. iBasri fc Batalhal [l990t iHartigan et a l. 1991) , with 
observations of accretors with low Av ([G ullbring et al.' 
2000.) ■ and with shock models (jCalv et & Gullbring 
19981) . 

Our extinction estimates are listed in Table 3. Any 

^ Obtained from |http://www-int. stsci.edu/ ~inr/ultracool.html| 
(Gizis et al. 2000ab, Kirkpatricic et al. 1999). 
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error in the Av for V927 Tau, MHO 7, CFHT 7, and 
2M0444-I-2512 will propogate to the remaining data. 
Our estimate s are consistent with prev i ous e xtinction 
estimates (Br iceno et al.l [20021 iLuhmaiJ |2004[ ). with a 
few exceptions. We estimate Av ^ 1.0 and ~ 0.7 
to 2M0455-I-3028 and 2M0439-H25 44, respec t ively com- 
pared with A J = 0.0 and 0.07 fromlLuhma^ (|2004f) . We 
measure Av — 2.0 ± 0. 7 to GM Tau, less than previous 
estimates of ^ 4.2 mag (jBriceno et al.ll2002l : iGuieu et al.l 



2M1207-3932 is in the TW Hya Association and iso- 
lated from molecular gas. We assume Av = to 2M1207- 
3932, and conflrm that Av < 1 by comparison to a tem- 
plate LRIS M8-dwarf standard. 

9.3. Stellar Radius and Mass 

Estimates for the stellar mass and radius are re- 
quired to calculate mass accretion rates. The spectral 
class is converted to t e mpera ture using the scale de- 
rived by iLuhman et al.l (|2003D . We then use our spec- 
tral types, extinction estimates, synthetic R-band pho- 
tometry (after correcting for the measured veiling), and 
Baraffe et al. (1998) evolutionary tracks to compute 
the stellar luminosity, mass, radius, and age (see Ta- 
ble 3). The calculated ages for the higher-mass sam- 
ple in Taurus are overestimated, which is consistent with 
known systematic proble ms in stellar evolutionary tracks 
(|Hillenbrand et al.|[2007[) 

A distance of 140 p c is assumed f or Ta urus stars 
(jKenvon et al.l 119941 ). ISlesnick et all (|2006[ ) identifled 
S0518-I-2327 as an accretor away from the main Tau- 
rus molecular cloud and classify it as an intermediate- 
age star based on the depth of Na I gravity indices. 
The distance to S0518-t-2327 is therefore more uncertain 
the most other Taurus stars. We adopt a parallax dis- 
tance of 54 ± 3 pc to 2M1207-3932 (Gizis et al. 2007; 
see also Biller fc Close 2007) and of 56 pc for TW Hya 
(jWichmann et al.lll998l ). 

Any binarity in our sample would reduce the ste llar lu- 
minosity and estimated radius. iKraus et al.l (|2006l ) found 
that GM Tau and MHO 5 are spatially unresolved to 
~ 0'.'05 and are either single stars or spectroscopic bina- 
ries. We assume our other targets are single stars with 
the exception of V927 Tau (see §9.4.2). 

9.4. Notes on Individual Objects 
9.4.1. 2M1207-3932 



2M1207-3932 is a M8-L dwarf binary ([ci^ [200l 
IChauvin et al.ll2004|) in the 10 Myr old T W Hya Associ- 
ation (jWebb et al. N1999I : lMamaiekll2005| ). The primary 
has an IR excess and ongoing accretion. The secondary is 
faint and undetected in our data. We observed 2M1207- 
3932 on both November 24 and February 8. The Febru- 
ary observation was deeper and obtained with the target 
at lower airmass. The flux from the November integra- 
tion is calibrated based on the strength of photospheric 
emission from 6000-8000 A measured in the February 
run. The blue and UV continuum emission was fainter 
in November than on February 8. Most of these param- 
eters have a higher confldence for 2M1207-3932 than the 
T aurus membe rs. 

IStelzer et all ()2007( ) applied the lNatta eFaLl (|2004D re- 
lationship for Ha 10% width to measure M ^ 10~^"'^ — 
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10-9-8 Mq yr-i from 2M1207-3932, which is two orders 
of magnitude larger than measured here. This differ- 
ence is not likely attributed to variability because the 
Ha equivalent width during our February observation is 
the largest yet detected from 2M1207-3932. During an 
epoch in 2003 with low M fjStelzer et al.ll2007[ ). both the 
equivalent width and Ha 10% width were small. The 
variable equivalent width in the Ca II A8542 line of < 0.3 
A (jStelzer et al.ll2007D corresponds to M < 2.1 x lO^i^ 
Mq yr~^ (Table 15), which confirms our low M. 

9.4.2. V927 Tau 

V927 Tau is a bi nary separated by ~ 0^2 7 
(iLeinert et al.l floM iWhite fc Ghe2 [200l . 
iMuzerolle etall (|2003f ) measured a broad 10% width 
of 290 km s~^ and only a small self-reversal in the 
symmetric Ha line profile from V927 Tau. They 
suggested that this Ha emission may be produced by 
chromospheric flaring since most Ha profiles produced 
by accretion are asymmetric. At the time no excess 
IR or sub-mm emission had been detected from cither 
component (Simon & Prato 1995; White & Ghez 2001i; 
[Andrews fc Williams! 120051) . However, iMcCabe et all 
(I2006D measured K- band magnitudes th at were the 
same as measured by I White fc GhezI (|2001[ ). but L-band 
emission that was brighter by 1 magnitude for both 
V927A and B. Such variability may suggest that the 
amount of warm dust in the disk is variable, and that 
perhaps accretion onto V927 Tau may sometimes be 
present. 

We detect some excess Balmer continuum and He I 
A5876 emission from V927 Tau relative to MHO 7. V927 
Tau and its template MHO 7 are slightly different spec- 
tral types, which may lead to some discrepancy in this 
comparison. Such emission could be interpreted as either 
chromsopheric or accretion-related. iHartigan fc KenvonI 
([2003) found weak excess Paschen continuum emission 
from the secondary, if the secondary is classified as an 
M3.5 star. V927 Tau remains ambiguous as to whether 
accretion is ongoing. 

9.4.3. 2M0i38+2611 

We obtain a short LRIS spectrum of 2M0438-I-2611, 
which iLuhman et al.l (|2007f ) describe as a young accret- 
ing star seen through an edge-on disk with a gray ex- 
tinction of 5.5 mag. The gray extinction cannot be mea- 
sured from the shape of the continuum emission, which 
is similar to that of 2M0444-I-2512 {Ay = mag). If 
we assume the two stars have the same radius, then we 
find Av = 2.8 mag to 2M0438-I-2611. That the opti- 
cal forbidden Unes are bright (see also iLuhmanI [200^ 
ILuhman et al.l [20071 ) indicates that some of the outflow 
emission is at least slightly extended and not absorbed 
by the edge-on disk, consistent with other accretors seen 
through edge on disks. 

We classify 2M0438-I-2611 as a moderate accretor 
based on the detection of optical forbidden lines and non- 
detection of the Ca II IR triplet. We use the H7 line 
luminosity and Table 12 to calculate Lace ~ 4 x 10~^ 
Lq, or M ^ 1.5 X 10^^^ Mq yr^^ (assuming same M 
and R are the same as 2M0444-I-2512). The upper limit 
for the Ca II A8542 equivalent width of < 1.2 A indi- 
cates Lace < 1-6 X lO^"* Lace- Wc do not use the M of 



2M0438-h2611 in any analysis. 

9.4.4. 2MO414+28II 

2M0414-I-2811 is about a magnitude br ighter and is 
bluer than UBV photometry obtained by iGrosso et al.l 
(12001 . which suggests a higher M during our observa- 
tion. 

9.4.5. MHO 5 

Emission in the [O I] A6300, 6363 and He I A5876 hues 
are detected with equivalent width typical of accretors. 
MHO 5 is therefore classified as an accretor, des pite a 
low Ha 10% width of 154 km s'^ (|Muzerolle et al.|[200a) . 
Ongoing accretion indicates that the weak excess Balmer 
continuum emission is likely produced by accretion, al- 
though we cannot rule out a chromospheric origin. 

9.4.6. CIDA 14 

The excess Balmer continuum emission onto CIDA 14 
is very weak, and only detected because the photospheric 
emission is very similar to that from the non-accretor 
MHO 7. We detect an equivalent width of 0.27 A in the 
[O I] line from CIDA 14, which is slightly lower than the 
0.4 A measured bv lBriceno et al.l (|1999f ) but larger than 
the < 0.1 A uppe r limit from high-resolution spectra 
by IMu zerolle e t al.l (|2003[ ). This low upper limi t may 
be un derestimated if the line is broad. iMuzerolle et al.l 
()2003() measured an Ha equivalent width of 289 km s"^. 
CIDA 14 is a likely accretor. 

9.4.7. KPNO 11 

Several He I lines that are undetected from the non- 
accretors are detected from KPNO 11 (Table 5). How- 
ever, the strength of the He I A5876 A line and the upper 
limit on the Na I D line luminosity are below the expected 
luminosity for the calculated M . These luminosities may 
suggest that the excess Balmer continuum emission and 
He I line s are chromospheric. The Ha 10% width of 182 
km s~^ (jMohantv et al.l I2005D is intermediate between 
accretors and non-accretors. 

9.4.8. GM Tau 

The flux ratio of Balmer line emission to Balmer con- 
tinuum emission from GM Tau is much smaller than mea- 
sured for any other source here, which requires a large 
optical depth in the slab. 

9.4.9. S0518+2327 

S0518-I-2327 is an accretor located away from the main 
Taurus Molecular Cloud (Slesnick et al. 2006). The dis- 
tance to this star is more uncertain the the other Taurus 
stars studied here. 

9.4.10. 2M0436+2351 

The U-band atmospheric correction for 2M0436-I-2351 
is somewhat uncertain. As a result, no slope for the 
Balmer continuum is listed in Table 5. The mass accre- 
tion rate for this star is also more uncertain than that of 
other stars in our sample. 
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9.4.11. DF Tau 



DF Tau is a O'.'l binary (jSchaefer et al.ll2006D that is 
unresolved in these spectra. When calculating the stellar 
radius we assumed that the luminosity of the two stars 
are equal. STIS observed DF Tau three times with the 
G430L grating, two of which are analyzed here. The 
third observations shows much weaker photospheric and 
continuum emission with an a ccretion rate 10 times lower 
than the other observations (jHerczeg et al.l [2006,1 . We 
do not analyze this third spectrum here because it was 
obtained with a narrower slit and could have included 



the secondary but not the primary. 

9.4.12. DG Tau 

We measure Af = 4.1 x 10^® Mq yr ~^, an order of 
magni tude less than that measured by iGullbring et al.l 
mm. The spect rum of DG Tau during the 
I Gullbring et al.l ()2000D observation was heavily-veiled 
and an order of magnitude brigther than that detected 
here. The different Af result from real variability. Strong 
emission is detected in forbidden lines and Fe I lines de- 
spite the lower M. 
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TABLE 1 
Observations with Keck //LRIS 



Dtar 


MJD° 


texp 


T T 

U 


b 


r> 


b 


V 


b 


T5t 

ti 


(7^ (mag) 


Abbreviation^ 


CIDA 1 


63, 


.318 


720 


18, 


.6 


18, 


.3 


17. 


2 


15.3 


0.6 


CIDA 1 


2MASS J04141188+2811535 


63, 


.329 


1200 


18. 


.2 


18, 


.2 


17. 


,6 


16.4 


0.6 


2M0414+2811 


MHO-Tau 5 


63, 


.382 


1620 


19, 


.2 


18, 


.3 


16. 


,9 


15.0 


0.4 


MHO 5 


MHO-Tau 6 


63, 


.405 


1260 


18, 


.5 


17, 


.8 


16. 


,6 


15.1 


0.4 


MHO 6 


GM Tau 


63, 


.422 


1050 


18, 


.7 


18, 


.7 


17. 


,8 


16.3 


0.4 


GM Tau 


2MASS J04390396+2544264 


63, 


.438 


2700 


22. 


.0 


21, 


.5 


20. 


1 


17.6 


0.4 


2M0439+2544 


2MASS J04442713+2512164 


63, 


.475 


2100 


21. 


.0 


20, 


.4 


19. 


,5 


17.3 


0.4 


2M0444+2512 


CIDA 14 


63, 


.498 


720 


18, 


.5 


17, 


.4 


15. 


,9 


14.4 


0.4 


CIDA 14 


SCH J0439016+2336030 


63, 


.508 


840 


18. 


,1 


17, 


.2 


15. 


,9 


14.6 


0.3 


S0439+2336 


V927 Tau 


63, 


.520 


150 


17, 


.3 


16, 


.1 


14. 


,8 


13.3 


0.3 


V927 Tau 


CFHT Tau 7 


63 


.525 


600 


20, 


.8 


20, 


.0 


18. 


,4 


16.3 


0.2 


CFHT 7 


MHO-Tau 7 


63, 


.532 


600 


19, 


.1 


18, 


.0 


16. 


,6 


15.0 


0.2 


MHO 7 


2MASS J04381486+2611399 


63, 


.540 


300 














20.1 


0.2 


2M0438+2611 


KPNO-Tau 11 


63, 


.552 


900 


20, 


,1 


19, 


.1 


17. 


,7 


15.9 


0.2 


KPNO 11 


2MASS J04414825+2534304 


63, 


.567 


2400 


22. 


.0 


21, 


.9 


21. 


1 


19.1 


0.2 


2M0441+2534 


2MASS J04554757+3028077 


63, 


.600 


1080 


19. 


.0 


17, 


.9 


16. 


,5 


14.9 


0.2 


2M0455+3028 


SCH J0518028+2327126 


63, 


.611 


1380 


18. 


.3 


18, 


.2 


17. 


,8 


16.6 


0.2 


S0518+2327 


2MASS J12073347-3932540'= 


63, 


.647 


840 


21, 


.1 


21, 


.3 


20. 


,5 


18.2 


0.2 


2M1207-3932 


2MASS J04362151+2351165 


139.218 


400 


20. 


.0 


19, 


.5 


18. 


,4 


16.9 


1.0 


2M0436+2351 


2MASS J12073347-3932540 


139.611 


1800 


20. 


.8 


20, 


.9 


20. 


,2 


18.0 


0.2 


2M1207-3932 



"540000+ 

''Synthetic photometry, calculated by convolving the spectra with Bessell UBVR filter transmission curves 
for the Nordic Optical Telescope at 'http : / / ww w. not .iac.es /instruments /filters/ 1 
'^Approximate flux uncertainty in magnitudes 
''Abbreviation for target used throughout this work. 
•^Used a 0'.'7 slit width. 



TABLE 2 

Archival HST/STIS G430L Observations 



Star 


MJD 




Slit 


PI 


Rcf. 


Abbreviation 


DF Tau 


51439.84 


144 


52 X 0.2 


Linsky 


2 


DF Tau (Sep. 99) 


TW Hya 


51671.88 


260 


52 X 0.2 


Linsky 


1 


TW Hya (May 00) 


TW Hya 


52474.58 


120 


52 X 0.2 


Johns-KruU 




TW Hya (Jul. 02) 


DF Tau 


52512.49 


120 


52 X 0.2 


Johns-KruU 




DF Tau (Aug. 02) 


RU Lupi 


51737.65 


120 


52 X 0.1 


Walter 


2 




V819 Tau"" 


51786.87 


360 


52 X 0.1 


Walter 


2 




DG Tau 


51839.52 


240 


52 X 0.1 


Walter 


2 




V836 Tau 


51954.36 


360 


52 X 0.1 


Walter 


2 




LkCa 7A'^ 


53039.71 


162 


52 X 0.2 


Johns-KruU 




LkCa 7 



1: Herczeg et al. (2004) 2: Herczeg et al. (2006) 

"We extract and use only the primary spectrum and do not discuss the faint secondary 
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TABLE 3 
Adopted Stellar Properties 



Star 


SpT 


log Age"' 


Av 


T 


L 


R 


M 


Ref. 






log "VT 


mag 


K 


Ln 


Rr, 


Mr, 




(^TT^ A 1 
OlUA 1 


IViO 


D.U 


o.U 




U.zi 


1.0 ( 


n 1 c; 
U. lo 


1 


\j iVi 1 clU 


IViO.O 


0.4 


9 n 


vSUy i 


U.U / U 


U.yo 


n 1 c; 
U. 10 


O 

Z 


01\/Tny1 1 A 1 OQ1 1 


IViO.O 


D.O 


1 n 
i.U 


zyoo 


U.U14 


n /I ft 
U.4D 


n nftc; 
U.Uoo 


Q 
O 


oUoioH-zoz 1 


1 Vi . 


^^ o 


i.U 


vSUy i 


n ni ft 
U.UID 


U.44 


n 1 n 
U.IU 


A 

4 


iVlxlC D 


IViO 


ft 1 
o. i 


1 n 
i.U 


Q 1 oc; 
vSizO 


U. 14 


1 Q 1 


n 1 7 
U. 1 / 







IViO 




n n 
U.U 


vSlZO 


U.Uoy 


1 no 
l.Uz 


n 1 7 


4 




IViO 




i.U 


oizo 


n ni Q 
U.Ulo 


n un 
U.oy 


nil 
U. 11 


ft 


Z iVl U 4 4 1 Z v5 4 


A/T7 




1 n 
i.U 


Z i OZ 


U.UUovj 


n 9c: 
U.zo 


n nQp; 
U.UoO 


Q 

O 


Z iVl U ^ ^"4: -p ^ 1 Z 


A/T7 9^^ 
IVi / . ZD 


D.D 


U.U 


Zooo 


n ni 9 
u.uiz 


u.^yiO 


u.uou 


o 


2M0439+2544 


M7.25 


6.5 


0.8 


2838 


0.015 


0.51 


0.050 


3 


2M1207-3932 


M8 


7.0 


0.0 


2550 


0.0022 


0.24 


0.024 


7 


MHO 5 


M6 


6.0* 


0.3 


2990 


0.098 


1.17 


0.10 


5 


2M0455+3028 


M4.75 


6.1 


1.0 


3161 


0.15 


1.31 


0.20 


3 


CIDA 14 


M5.5 


6.0 


0.3 


3058 


0.16 


1.44 


0.12 


1 


KPNO 11 


M5.75 


6.5 


0.3 


3023 


0.037 


0.70 


0.10 


8 


V927 Tau'' 


M5 


6.0* 


0.4 


3125 


0.22 


1.60 


0.15 


9 


MHO 7 


M5.5 


6.3 


0.1 


3058 


0.076 


0.99 


0.14 


5 


CFHT 7 


M6.5 


6.5 


0.0 


2935 


0.025 


0.61 


0.07 


10 


HST/STIS Sample 


TW Hya 


K7 


7.6 


0.0 


4060 


0.17 


0.83 


0.77 


11 


RU Lup 


K7 


7.0 


0.1 


4060 


0.42 


1.3 


1.0 


12 


DF Tau* 


M2 


6.0* 


0.6 


3560 


1.45 


3.2 


0.66 


13 


DG Tau 


K6 


7.5 


1.3 


4205 


0.28 


1.0 


0.85 


13 


V836 Tau 


K7 


7.1 


1.1 


4060 


0.32 


1.15 


0.95 


14 



Bricenoetal 

^^gren 

Luhman (200^ 
Slesnick et al. (2006) 
Briceno et al. (1998) 
Luhman et al. (2006) 





Gizis (2002) 

''I Myr is the youngest age of ^araffe_e^_aL (_1998) models 

Ages marked with a "*" are upper limits. 
''Binary, assumes both components have same M, L, and T. 



GuieuetaL 
Ru cinski fc Krautter (1983) 
Henize (1954) 
Smak (1964) 
Feigelson fc Kriss (1981) 
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TABLE 4 

Veiling Measurements of M-dwarfs from Keck spectra" 



Star 


Template 


4227 


4955 


5450 


6190 


6880 


7100 


7800 


8600 


CIDA 1 


V927 Tau 


> 2 


0.98 


0.75 


0.94 


0.42 


0.63 


0.38 


0.56 






0.85 


1.4 


2.1 


3.6 


6.1 


8.0 


10.9 


25.0 


GM Tau 


MHO 7 


> 2 


1.84 


1.55 


3.63 


0.88 


0.79 


0.27 


0.26 






0.8 


1.2 


1.6 


2.7 


3.7 


3.5 


2.9 


4.7 


2M0414+2811 


CFHT 7 


> 2 


> 3 


> 3 


> 3 


1.52 


1.59 


0.76 


0.16 






1.6 


2.3 


2.5 


3.0 


4.2 


4.3 


4.7 


2.9 


S0518+2327 


MHO 7 


> 2 


1.62 


0.79 


0.87 


0.38 


0.38 


0.14 


0.40 






1.6 


1.7 


1.4 


1.5 


1.8 


1.9 


1.5 


5.0 


MHO 6 


V927 Tau 


0.49 


0.05 


0.03 


0.07 


0.02 


0.01 


-0.05 


0.02 






0.68 


< 0.41 


< 0.56 


< 0.55 


< 0.75 


< 1.6 


< 1.7 


< 2.6 


sni.'in-2.'!:iri 


Vn27 Tail 


0.18 


0.0.") 


0.0.1 


n.02 


0.02 


0.01 


0.04 


0.1. 'i 








^ 0.7(i 


< 1.0 


^ 1.0 


^ 1.2 


^ 2.4 


< 2.() 


^ ;i.8 


2M0436+2351 


V927 Tau 


0.36 


0.08 


0.04 


-0.00 


-0.02 


0.02 


-0.07 


-0.12 






0.100 


0.135 


0.092 


-0.001 


-0.120 


0.134 


-0.820 


-1.817 


2M0455+3028 


V927 Tau 


0.04 


0.02 


0.03 


0.04 


0.02 


0.02 


-0.01 


-0.03 






< 0.095 


< 0.43 


< 0.60 


< 0.66 


< 0.88 


< 1.9 


< 2.2 


< 3.7 


2M0441+2534 


2M0439+2544 


> 2" 


> 3 


2.4 


4.1 


0.50 


0.61 


0.29 


-0.01 






0.04 


0.06 


0.084 


0.14 


0.17 


0.20 


0.24 


< 0.085 


2M0444+2512 


2M0439+2544 


2.06'' 


0.11 


0.01 


0.06 


0.04 


0.05 


-0.03 


0.03 






0.079 


0.060 


< 0.04 


< 0.04 


< 0.07 


< 0.17 


< 0.21 


< 0.45 


2M1207-3932 


2M0439+2544 


0.29'' 


0.32 


0.04 


0.59 


0.10 


0.07 


-0.12 


-0.01 


(Feb.) 




0.018 


0.069 


< 0.017 


< 0.015 


< 0.038 


< 0.084 


< 0.11 


< 0.28 


2M1207-3932 


2M0439+2544 


- 


0.13 


0.05 


0.19 


0.03 


0.01 


-0.10 


0.01 


(Nov.) 




— 


0.028 


< 0.015 


< 0.011 


< 0.032 


< 0.080 


< 0.13 


< 0.24 


MHO 5 


CFHT 7 


-0.01 


-0.02 


0.01 


0.08 


0.03 


-0.01 


-0.10 


-0.03 






< 0.06 


< 0.34 


< 0.42 


< 0.41 


< 0.76 


< 1.7 


< 2.0 


< 3.6 


CIDA 14 


MHO 7 


-0.02 


-0.00 


-0.00 


0.01 


0.01 


-0.01 


-0.04 


0.05 






< 0.16 


< 0.80 


< 1.0 


< 1.0 


< 1.4 


< 3.1 


< 3.4 


< 5.6 


KPNO 11 


MHO 7 


0.06 


0.00 


-0.06 


-0.04 


-0.02 


-0.03 


-0.05 


-0.03 






< 0.027 


< 0.16 


< 0.20 


< 0.21 


< 0.32 


< 0.74 


< 0.89 


< 1.6 


2M0438+2611 


2M0439+2544 




-0.04 


-0.29 


0.14 


0.12 


0.04 


0.03 


0.01 








< 002 


< II 0112 


< II Oil 1 


< II OIK) 


< 01-4 


< 018 


< 0.36 








l'li()(o,-,| 


)lic;rif UViii] 


jial.rs 










2M0439+2544 


CFHT 7" 


0.54 


0.06 


0.03 


0.13 


-0.02 


-0.02 


-0.08 


-0.10 






0.011 


< 0.013 


< 0.020 


< 0.023 


< 0.052 


< 0.13 


< 0.17 


< 0.38 


V927 Tau 


MHO 7° 


0.04 


0.14 


0.10 


0.06 


0.07 


0.09 


0.10 


0.36 






< 0.54 


< 2.3 


< 2.9 


< 3.1 


< 3.9 


< 8.2 


< 8.8 


< 13.3 


CFHT 7 


MHO 7= 


0.09 


-0.07 


-0.09 


-0.07 


-0.07 


-0.07 


-0.04 


-0.22 






< 0.009 


< 0.075 


< 0.10 


< 0.11 


< 0.20 


< 0.48 


< 0.64 


< 1.3 


MHO 7 


V927 Tau= 


-0.08 


0.08 


0.10 


0.08 


0.08 


0.07 


-0.09 


-0.27 






< 0.093 


< 0.44 


< 0.55 


< 0.54 


< 0.77 


< 1.7 


< 1.9 


< 3.4 



'^Veiling (accretion continuum fiux divided by photosphcric continuum flux) at listed wavelength is given on top line 
accretion continuum flux (10"^*^ erg cm~^ s""'^ ^~^) given on bottom line for each source. 
Flux upper limits are marked with a < and are always calculated for rx < 0.05. 
Stars listed in order of M (see Table 11). 

^'CFHT 7 is used as a templalc instead of 2M0439+2544 at 4227 A. 

■^Mismatched spectral template. 



TABLE 5 

Veiling Measurements for higher-mass accretors observed with HST^ 



Star 


Template 


4000 


4600 


5400 


DF Tau 


LkCa 7 


> 5 


1.8 


1.3 


(Sep. 99) 




310 


310 


275 


DF Tau 


LkCa 7 


> 5 


0.73 


0.6 


(Aug. 02) 




150 


126 


114 


RU Lup 


V819 Tau 


> 5 


1.7 


2.3 






840 


590 


710 


DC Tau 


V819 Tau 


3.0 


1.4 


1.7 






37 


59 


64 


TW Hya 


V819 Tau 


2.2 


0.58 


0.36 


(May 00) 




510 


450 


330 


TW Hya 


V819 Tau 


1.4 


0.69 


0.72 


(Jul. 02) 




360 


460 


470 


V836 Tau 


V819 Tau 


0.23 


< 0.1 


< 0.08 






3.1 


< 4.9 


< 4.0 



'^Format same as in Table 4 



Accretion onto very low-mass stars and brown dwarfs 



TABLE 6 

Observed properties of the accretion continuum 



Star 


ja 


J*" 




^4966/^4^227 


Accretors in Keck I sample 


CIDA 1 


3.4 


5.6 ±0.7 


4.1 ±2.0 


0.85 


GM Tau 


3.5 


5.9 ±0.7*^ 


6.1 ±1.5 


0.9 


2M0414+2811 


2.7 


3.10 ±0.17 


8.1 ±1.0 


1.14 


S0518+2327 


2.1 


2.35 ±0.15 


2.4 ±2.0 


0.83 


MHO 6 


1.10 


3.9 ±0.8 


7.7 ±1.5 


< 0.5 


S0439+2336 


0.76 


3.1 ± 1.5 


6.2 ±1.1 


< 1.2 


2M0436+2351 


1.48 


3.8 ± 1.2 






2M0441+2534 


2.8 


3.75 ± 0.3 


7.2 ± 1.0 


1.2 


2MU444+2512 


1.24 


3.1 ±0.7 


6.6 ± 1.2 


0.76 


2M0439+2544 


1.37 


> 3 


16 ±6 




2M1207-3932 


5.9 


15 ±4 


8.7 ± 1.0 




2M1207-3932 


3.4 


7±2 


10.7 ±2.0 




Possible 


accretors or non-accretors in Keck I sample 


2M0455+3028 


0.46 


> 1.3 






MHO 5 


0.46 


> 1.5 






CIDA 14 


0.35 








KPNO 11 


0.41 








V927 Tau 


0.38 








MHO 7 


0.34 








CFHT 7 


0.36 








Accretors in HST sample 


DF Tau (Sep. 99) 


1.58 


1.9 ±0.1 


5.4 ±0.4 




DF Tau (Aug. 02) 


1.95 


2.8 ±0.3 


5.0 ±0.4 




RU Lup 


1.72 


1.9 ±0.2 


2.9 ±0.5 




DG Tau 


1.35 


2.0 ±0.1 


5.2 ±0.6 




TW Hya (May 00) 


1.66 


2.3 ±0.2 


4.3 ±0.6 




TW Hya (Jul. 02) 


2.0 


3.1 ±0.5 


4.8 ±0.4 




V836 Tau 


0.60 










Non- 


accretors in 


HST sample 




LkCa 7 


0.44 








V819 Tau 


0.46 









"Observed flux ratio, ^4000 7^3600 

''Flux ratio, i^4000 /^3600 , after photospheric subtraction and dereddening 
■^Slope of Balmer continuum 

''Flux ratio as proxy for slope of Paschen continuum. 
''Measured at alternate wavelength due to strong line blending 
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TABLE 7 

Balmer Line Equivalent Widths and Fluxes 



Star 


Hq 


H/3 


H7 


H(5 


H6 


H7 


H8 


H9 


HIO 


Hll 


H12 


CIDA 1 


232 


134 


80 


64 


41 


49 


37 


27 


18 


17 


10 




31.2 


3.0 


1.08 


0.65 


0.38 


0.38 


0.29 


0.24 


0.21 


0.20 


0.14 


GM Tau 


281 


68 


39 


30 


8.8 


26 


13 


11 


12 


11 


9.1 




14.6 


1.11 


0.41 


0.24 


0.12 


0.21 


0.13 


0.11 


0.15 


0.14 


0.12 


2M0414+2811 


280 


103 


72 


58 


30 


32 


23 


17 


10.6 


8.7 


5.6 




15 


2.5 


1.3 


0.97 


0.57 


0.53 


0.39 


0.32 


0.26 


0.22 


0.15 


S0518-L2.'527 


r-,a 


f)7 


60 


rii 


28 


27 


19 


9,8 


5.3 


4,7 


1.9 




2.') 


l.Til 


1.27 


1.11 


U.li!) 


0.()2 


().4:i 


0.2(i 


O.IS 


O.Ki 


0,069 


MHO 6 


55 


31 


37 


28 


27 


27 


21 


16 


12.1 


9.6 


5.3 




10.2 


1.6 


1.00 


0.73 


0.45 


0.41 


0.32 


0.26 


0.19 


0.16 


0.097 


S0439+2336 


13 


18 


24 


18 


21 


16 


14 


11 


8.9 


6.2 


2.7 




3.9 


1.55 


1.24 


0.90 


0.60 


0.47 


0.39 


0.31 


0.22 


0.16 


0.078 


2M0436+2351 


28 


47 


74 


62 


46 


44 


32 


16 


10.2 


9.2 


5.7 




0.97 


0.48 


0.38 


0.29 


0.19 


0.16 


0.106 


0.064 


0.046 


0.046 


0.028 


2M0441+2534 


586 


228 


158 


135 


72 


70 


47 


29 


14.6 


9.1 


3.9 




1.7 


0.16 


0.085 


0.061 


0.036 


0.029 


0.020 


0.014 


0.0101 


0.0070 


0.0035 


2M0444+2512 


100 


126 


94 


62 


49 


27 


20 


18 


10.3 


B 


B 




1.7 


0.34 


0.15 


0.085 


0.058 


0.031 


0.022 


0.018 


0.012 


- 


- 


2M0439+2544 


158 


93 


93 


87 


73 


50 


44 


36 


36 


18 


6.4 




1.7 


0.12 


0.054 


0.035 


0.022 


0.016 


0.012 


0.0093 


0.0079 


0.0056 


0.0027 


2M1207-3932 


396 


341 


286 


323 


142 


146 


132 


89 


55 


27 


16 


(Feb.) 


3.7 


0.53 


0.27 


0.20 


0.124 


0.099 


0.084 


0.066 


0.050 


0.037 


0.027 


2M1207-3932 


200 


210 


234 


132 


- 


- 


- 


- 


- 


- 


- 


(Nov.) 


1.21 


0.19 


0.119 


0.088 


0.071 


0.052 


0.042 


0.037 


0.024 


0.0106 


0.008 


2M0438+2611 


99 


158 


622 


- 


- 


- 


- 


- 


- 


- 


- 




0.16 


0.023 


0.013 


0.0056 


0.0052 


< 0.0036 


< 0.0036 


< 0.0024 


< 0.0014 


< 0.0014 


< 0.0014 


2M0455+3028 


9.6 


9.2 


8.6 


6.1 


7.5 


7.1 


4.0 


4.5 


3.0 


2.0 


< 1.0 




2.12 


0.45 


0.26 


0.15 


0.095 


0.082 


0.051 


0.054 


0.029 


0.018 


< 0.012 


MHO 5 


43 


30 


19 


13.4 


13.1 


10.4 


6.5 


5.7 


5.2 


2.8 


2.1 




6.6 


0.83 


0.34 


0.20 


0.119 


0.084 


0.058 


0.049 


0.029 


0.020 


0.015 


CIDA 14 


10.2 


7.5 


3.9 


3.4 


3.8 


4.5 


1.4 


2.2 


< 5.4 


< 3.8 


< 3.8 




3.3 


0.54 


0.21 


0.15 


0.083 


0.089 


0.032 


0.048 


0.07 


< 0.05 


< 0.05 


KPNO 11 


13.6 


12.5 


9.5 


7.6 


8.4 


7.4 


3.3 


3.6 


3.5 


< 1.6 


< 1.6 




1.03 


0.18 


0.087 


0.054 


0.033 


0.028 


0.014 


0.014 


0.010 


< 0.005 


< 0.005 


V927 Tau 


10.7 


7.9 


3.7 


3.2 


3.6 


3.9 


1.6 


2.2 


2.1 


< 1.7 


< 1.7 




10.2 


1.9 


0.71 


0.46 


0.26 


0.25 


0.12 


0.15 


0.11 


< 0.08 


< 0.08 


MHO 7 


7.4 


5.4 


2.5 


1.9 


2.0 


< 2.5 


< 1.5 


1.6 


< 5.7 


< 1.7 


< 2.0 




1.31 


0.22 


0.079 


0.049 


0.026 


< 0.03 


< 0.02 


0.021 


< 0.04 


< 0.012 


< 0.015 


CFHT 7 


15 


14 


14 


8.5 


7.4 


5.3 


2.5 


3.5 


6.0 








0.63 


0.084 


0.042 


0.021 


0.0115 


0.0082 


0.004 


0.005 


0.0005 


< 0.003 


< 0.003 



Table 7 (continued): HST/STIS Sample 



Star Ha H/3 


H7 


H5 


H6 


H7 


H8 


H9 


HIO 


Hll 


H12 


DP Tau 


29 


27 


26 


22 


17 


13 


8.4 


5.2 


3.4 


2.9 


(Sep. 99) 


150 


109 


90 


68 


54 


42 


30 


19 


13 


11.5 


DF Tau 


26 


20 


17.7 


16 


15 


10.9 


7.8 


6.1 


3.9 


4.2 


(Aug. 02) - 


80 


44 


32 


24 


24 


16 


13.6 


11.3 


7.6 


7.9 


RU Lupi 


37 


16 


14.5 


7.3 


8.7 


6.6 


9.3 


4.6B 


4.4 


3.0 




360 


160 


113 


68 


74 


57 


80 


51B 


48 


33 


DG Tau 


24 


14 


11.6 


8.8 


11.8 


6.9 


4.2 


3.6 


2.5 


5.3 




29 


11.5 


7.0 


4.2 


5,1 


3.1 


2.1 


1.9 


1.4 


2.5 


TW Hya 


38 


28 


28 


24 


20 


14 


9.5 


7.0 


4.1 


5.0 


(May 00) 


470 


280 


230 


160 


137 


104 


76 


57 


37 


42 


TW Hya 


77 


48 


42 


33 


26 


18 


12 


7.6 


5.3 


3.6 


(Jul. 02) 


880 


420 


290 


190 


150 


113 


84 


58 


43 


30 


V836 Tau 


17 


8.8 


9.4 


13 


9.6 


7.5 


1.0 


6.2 


5.1 


3.4 




9.6 


3.1 


2.0 


1.3 


0.84 


0.70 


0.85 


0.45 


0.35 


0.25 


V819 Tau 


< 0.8 


< 2 


< 1.8 


< 4 


< 7 


< 2 


< 6 


< 6 


< 6 






< 1.3 


< 1.6 


< 0.9 


< O.i 


i < 0.6 


< 0.4 


< 1.0 


< 0.6 


< 0.8 




LkCa 7 


1.8 


< 3 


2.0 


3.2 


4.7 


1.9 


< 5.5* 


< 5 


< 7 






4.8 


< 6 


2.2 


1.5 


1.6 


0.86 


< 2.2* 


< 1.5 


< 2.5 




'^Equivalent width (-A) in 


first row 


flux 


(10-" 


erg 


s-^) in 


second row for each star. 





B=blend, *=inarginal detection, ID=questionable line identification 



Accretion onto very low-mass stars and brown dwarfs 



TABLE 8 

He Line Equivalent Widths and Fluxes" 



Star 


A4026 


A4470 


A4922 


A5016 


A5876 


A6678 


A7065 


He II A4686 


Keck //LRIS sample of low-mass stars 


CIDA 1 


4.3B 


5.3 


8.3 


10.4 


9.6 


1.9 


2.0 


1.3 




0.035B 


0.086 


0.22 


0.26 


0.49 


0.22 


0.31 


0.23 


GM Tau 


2.5 


5.0 


13.9 


16.1 


4.6 


1.1 


0.5 


0.2 




0.019B 


0.055ID 


0.23 


0.26 


0.12 


0.052 


0.030 


0.015 


2M0414+2811 


2.8 B 


4.2 


8.2 


9.6 


7.3 


2.4 


1.5 


1.1 




0.038 B 


0.085 


0.20 


0.24 


0.23 


0.11 


0.084 


0.062 


S0518+2327 


3.5 


5.7 


2.8 


3.4 


7.7 


2.4 


0.9 


1.5 




0.064 


0.12 


0.071 


0.078 


0.19 


0.087 


0.041 


0.033 


MHO 6 


1.9 


2.1 


0.6 


0.9 


2.9 


0.5 


0.6 


0.5 




0.047 


0.082 


0.041 


0.044 


0.19 


0.060 


0.10 


0.023 


S0439+2336 


1.2 


1.2 


0.3 


0.4 


1.5 


0.4 


< 0.2 


0.1 




0.054 


0.094 


0.033 


0.038 


0.18 


0.084 


< 0.07 


0.011 


2M0436+2351 


3.8 


3.8 


0.95 


< 0.75* 


2.9 


1.7 


2.3 


0.84 




0.018 


0.028 


0.012 


< 0.009* 


0.052 


0.036 


0.103 


0.0076 


2M0441+2534 


8.1 


16.2 


6.8 


7.5 


29 


8.5 


3.6 


2.1 




0.0030 


0.0084 


0.0051 


0.0058 


0.041 


0.021 


0.013 


0.0012 


2M0444+2512 


3.9 


7.2 


1.5 


4.9 


12.4 


1.4 


0.7 


2.1 




0.0053 


0.014 


0.0063 


0.014 


0.054 


0.014 


0.013 


0.0049 


2M0439+2544 


4.6 


8.7 


1.0 


2.3 


11.2 


1.8 


0.5 


1.5 




0.0022 


0.0062 


0.0022 


0.0031 


0.028 


0.013 


0.0075 


0.0015 


2M1207-3932 


8.8 


22 


2.5 


3.6 


33 


4.2 


1.6 


3.1 


(Feb.) 


0.0074 


0.021 


0.0056 


0.0057 


0.075 


0.022 


0.014 


0.0039 


2M1207-3932 


< 6.5 


267 


3.7 


4.5 


21 


3.5 


1.2 


< 3.1 


(Nov.) 


< 0.004* 


0.017 


0.0069 


0.0053 


0.038 


0.015 


0.010 


< 0.0034 


2M0438+2611 










< 10 


< 2 


< 1.8 






< 0.003 


< 0.003 


< 0.002 


< 0.002 


< 0.004 


< 0.0015 


< 0.0025* 


< 0.002 


2M0455+3028 


0.3 


< 0.17 


0.1 


0.1 


0.4 


< 0.13 


< 0.3 


< 0.08 




0.0081 


< 0.007* 


0.0070 


0.0056 


0.030 


< 0.02* 


< 0.06 


< 0.004 


MHO 5 


< 0.8* 


< 1.0* 


< 0.4* 


< 0.4* 


1.3 


< 0.25* 


< 0.16* 


< 0.25 




< 0.014 


< 0.020 


< 0.022 


< 0.042 


0.0555 


< 0.025 


< 0.025 


< 0.006 


CIDA 14 


< 0.3 


< 0.18 


< 0.1 


n.i 


n.3 


< 0.07 


< 0.1 


< o.m 




o.ui:) 


< 0.012 


< 0.012 


< O.OUi) 


u.u:i7 


< O.Oll) 


< 0.04 


< 0.007 


KPNO 11 


0.4 


0.4 


0.2 


< 0.18 


0.4 


0.1 


< 0.2 


< 0.10 




0.0033 


0.0053 


0.0059 


< 0.003 


0.0092 


0.0061 


< 0.015 


< 0.0015 


V927 Tau 


< 0.15 


< 0.4 


< 0.15 


< 0.02 


0.30 


< 0.2 


< 0.2 


< 0.25 




< 0.02 


< 0.08 


< 0.06 


< 0.06 


0.11 


< 0.12 


< 0.18 


< 0.06 


MHO 7 


< 0.18 


< 0.33 


< 0.25 


< 0.30 


< 0.20 


< 0.27 


< 0.17 


< 0.2 




< 0.005 


< 0.013 


< 0.02 


< 0.014 


< 0.012 


< 0.035 


< 0.035 


< 0.008 


CFHT 7 


< 0.45 


< 0.4 


0.20 


< 0.26 


0.35 


0.13 


< 0.3 


< 0.4 




< 0.0015 


< 0.002 


0.0027 


< 0.002 


0.0042 


0.0035 


< 0.015 


< 0.002* 



Table 8 (conlinuod) 



Sl.ar 


Ain2() 


A 1171) 


A 1922 


A5()l() 


A.lN7(i 


A(i()7,s 


A70(i.-) 


Ik! 11 A4()86 


HST/STIS sample of higher mass stars 


DF Tau 


1.9 


2.8 


2.0 


1.3 








0.7 


Sep. 99 


6.1 


12.4 


10.8 


6.6 








3.3 


DF Tau 


1.7 


1.3 


2.4 


2.6 








0.4 


Aug. 02 


2.8 


3.5 


8.3 


7.9 








1.13 


RU Lupi 


3.3B 


3.3 


9.5 


9.5 








< 0.35* 




25B 


33 


90 


92 








< 3.2* 


DG Tau 


1.3 


1.7 


3.2 


4.4 








< 0.5 




0.65 


1.5 


3.9 


5.0 








< 0.5 


TW Hya 


2.2 


2.8 


0.9 


1.0 








1.1 


May 00 


16.8 


29.8 


10.9 


11.8 








14.1 


TW Hya 


2.3 


2.9 


1.5 


1.2 








1.0 


Jul. 02 


14.9 


28 


17 


12.8 








11.4 


V836 Tau 


1.1 


0.6 


0.7 


0.3 








0.5 




0.18 


0.26 


0.40 


0.19 








0.27 


V819 Tau 


< 2.5 


< 1.5 


< 1.4 


< 0.7 








< 1.0 




< 1.0 


< 1.5 


< 2.0 


< 1.0 








< 1.5 


LkCa 7 


< 2.8 


< l.i 


< 1.2 


< 0.4 








< 1.2 




< 2,(1 


. 2 . 


< :i,n 


l.n 








< :!,() 



'^Equivalent width and line flux format as in Table 4. 

B=blend, *=marginal detection, ID=questionable line identification 
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TABLE 9 

Ca and Na Line Equivalent Widths and Fluxes 



Star 


Ca II K 


Ca II H"^ 


Ca II A8498 


Ca II A8542 


Ca II A8662 


Na I A5893 


O I A7773 


O I A8446 


CIDA 1 


175 


160 


14.7 


18.9 


14.4 


10.4 


1.7 


2.5 




1.49 


1.32 


7.8 


10.4 


9.1 


0.48 


0.54 


1.36 


GM Tau 


106 


135 


28 


29 


21 


15.4 


1.14 


3.4 




1.35 


1.15 


4.4 


4.5 


4.1 


0.41 


0.114 


0.55 


2M0414+2811 


116 


142 


39 


40 


30 


11 


1.24 


2.0 




1.9 


2.3 


4.9 


5.4 


5.1 


0.35 


0.108 


0.28 


S0518-L2.'?27 


20 


r-<9 


1.73 


2.1 


1.04 


4.3.1 


2.30 


1,23 




(1.4:! 


1.21 


(1.22 


(1.27 


(i.;i2 


(1.(11):! 


0.1!) 


0.15 


MHO 6 


30 


43 


< 0.9 


< 1.5 


< 0.4 


0.98 


< 0.4 


< 0.5 




0.45 


0.92 


< 0.4 


< 0.7 


< 0.3 


0.048 


< 0.14 


< 0.24 


S0439+2336 


23 


35 


< 0.08 


< 0.15 


< 0.09 


< 0.36 


0.16 


< 0.17 




0.59 


1.21 


< 0.06 


< 0.11 


< 0.08 


< 0.04 


0.083 


< 0.13 


2M0436+2351 


41 


79 


1.7 


1.6 


1.1 


2.2 


1.1 


< 0.5 




0.13 


0.35 


0.15 


0.16 


0.12 


0.039 


0.065 


< 0.5 


2M0441+2534 


53 


160 


6.6 


6.4 


5.0 


11.4 


1.14 


1.29 




0.021 


0.062 


0.094 


0.107 


0.123 


0.016 


0.0090 


0.019 


2M0444+2512 


49 


83 


< 0.6* 


0.45 


0.29 


20 


< 0.4 


< 0.6 




0.053 


0.101 


< 0.04* 


0.030 


0.028 


0.068 


< 0.016 


< 0.05 


2M0439+2544 


38 


93 


< 0.8 


< 1.1 


< 0.8 


4.8 


< 0.4 


< 1.1 




0.013 


0.037 


< 0.05 


< 0.07 


< 0.08 


0.0088 


0.013 


< 0.08 


2M1207-3932 


36.36 


196.32 


< 0.7 


< 1.6 


< 0.8 


5.4 


2.02 


1.95 


(Feb.) 


0.024 


0.178 


< 0.03 


< 0.076 


< 0.05 


0.013 


0.048 


0.082 


2M1207-3932 


90 


200 


< 1.3 


< 0.9 


< 0.9 


4.7 


< 0.9 


< 1.3 


(Nov.) 


0.012 


0.071 


< 0.06 


< 0.05 


< 0.03 


0.008 


< 0.02 


< 0.06 


2M0438+2611 


- 


34 


< 1.3 


< 1.2 


< 0.8 


< 17 


< 0.9 


< 0.9 




- 


0.0053 


< 0.008 


< 0.008 


< 0.007 


< 0.007 


< 0.003 


< 0.005 


2M0455+3028 


30 


26 


< 0.3 


< 0.3 


0.31 


0.16 


< 0.11 


< 0.35 




0.37 


0.43 


< 0.21 


< 0.23 


0.28 


0.009 


< 0.049 


< 0.25 


MHO 5 


48 


42 


< 1.3 


< 0.8 


< 0.7 


1.6 


< 0.2 


< 0.5 




0.43 


0.51 


< 0.8 


< 0.5 


< 0.6 


0.046 


< 0.07 


< 0.3 


CIDA 14 


14.3 


11.9 


< 0.10 


< 0.15 


< 0.10 


0.86 


< 0.12 


< 0.15 




0.30 


0.34 


< 0.10 


< 0.14 


< 0.13 


0.057 


0.080 


< 0.22 


KPNO 11 


32 


26 


< 0.25 


< 0.31 


0.28 


< 0.20 


< 0.12 


< 0.26 




0.125 


0.143 


< 0.07 


< 0.09 


0.123 


< 0.004 


< 0.021 


< 0.08 


V927 Tau 


24 


18.4 


< 0.21 


< 0.17 


< 0.13 


< 0.20 


< 0.13 


< 0.34 




1.55 


1.63 


< 0.5 


< 0.4 


< 0.40 


< 0.06 


< 0.31 


< 0.9 


MHO 7 


13 


10 


< 0.2 


< 0.2 


< 0.2 


< 0.2 


< 0.2 


< 0.4 




0.17 


0.18 


< 0.10 


< 0.11 


< 0.14 


< 0.011 


< 0.08 


< 0.27 


CFHT 7 


27 


21 


< 0.4 


< 0.4 


0.21 


< 0.4 


< 0.2 


< 0.6 




0.043 


0.047 


< 0.09 


< 0.09 


0.064 


< 0.004 


< 0.025 


< 0.14 



■^Equivalent width and line flux format as in Table 4. 

B— blend, *— marginal detection, ID=questionable line identification 

"^Blended with He 



Accretion onto very low-mass stars and brown dwarfs 



TABLE 10 

Wind Line Equivalent Widths and Fluxes 



Star 


O I A6300 


O I A6363 


O I A5577 


S II A6716 


S II A6731 


Fo II A7155 


S II A4069 


S II A4076 


CIDA 1 


1.17 


<0.32 


0.48 


<0.28 


0.21 


<0.30* 


<3.5B 


<2.7B 




0.084 


<0.029 


0.0096 


<0.030 


0.023 


<0.042* 


<0.034B 


<0.028B 


GM Tau 


5.4 


2.0 


1.9 


0.34 


0.64 


<0.27* 


<7.5B 


<5.0B 




0.19 


0.077 


0.052 


0.015 


0.030 


<0.015* 


<0.069B 


<0.040B 


2M0414+2811 


3.2 


1.17 


1.01 


0.42 


0.17 


0.15 


<2.8B 


<2.8B 




0.12 


0.047 


0.033 


0.020 


0.0078 


0.0081 


<0.047B 


<0.049B 


S0518-i-2.'?27 


1.02 


0.44B 


1.30 


<0,21* 


<n.i7 


0.25 


<0.63 


<1.47B 




O.U2tS 


U.U14B 


().():!:> 


<,().OU()IJ* 


<u.uo:)7 


O.0U9IJ 


<0.()127 


<^U.U:!UB 


MHO 6 


0.88 


0.17 


0.18 


<0.12 


<0.31 


<0.24 


6.0 


2.2 




0.075 


0.021 


0.021 


<0.014 


<0.036 


<0.032 


0.15 


0.051 


S0439+2336 


0.46 


<0.18 


0.17 


<0.21 


<0.15 


<0.18 


0.68 


0.53 




0.068 


<0.037 


0.028 


<0.040 


<0.028 


<0.041 


0.031 


0.024 


2M0436+2351 


< 1.0 


< 0.3 


< 1.4 


< 1.5 


< 0.2 


< 0.7 


< 2.0 


< 2.0 




< 0.015 




< 0.019 


< 0.03 


< 0.005 


< 0.02 


< 0.011 


< 0.011 


2M0441+2534 


8.3 


2.3 


4.3 


<1.02* 


<0.67* 


<0.26 


<5.29 


<2.60 




0.015 


0.0048 


0.0011 


<0.0026* 


<0.0017* 


<0.0008 


<0.0023 


<0.0012 


2M0444+2512 


62 


13.2 


6.6 


16 


29 


0.68 


413 


129 




0.40 


0.13 


0.038 


0.13 


0.27 


0.0094 


0.56 


0.16 


2M0439+2544 


1.51 


<0.36 


1.75 


<0.93 


<0.18 


<0.18 


<1.39 


<1.66 




0.0060 


<0.0023 


0.0071 


<0.0059 


<0.0012 


<0.0016 


<0.0006 


<0.0008 


2M1207-3932 


< 2.3 


< 0.7 


2.4 


< 0.47 


< 0.44 


< 0.38 


< 3.30 


7.1 


(Feb.) 


< 0.006 


< 0.003 


0.0059 


< 0.0022 


< 0.0022 


< 0.0024 


< 0.0026 


0.0062 


2M1207-3932 


< 5.0 


<0.91 


1.6 


<0.53 


<0.36 


<0.36 


<8.14 


<7.46 


(Nov.) 


< 0.013 


<0.0035 


0.0015 


<0.0021 


<0.0015 


<0.0020 


<0.0047 


< 0.0048 


2M0438+2611 


98 


20 


20 


5.0 


8.6 


<2.3 


- 


- 




0.048 


0.0149 


0.0052 


0.0037 


0.0072 


<0.0023 


0.0106 


0.0018 


2M0455+3028 


<0.06 


<0.09 


0.27 


<0.20 


<0.11 


<0.12 


<0.22 


<0.14 




<0.0066 


<0.0124 


0.0112 


<0.028 


<0.0147 


<0.021 


<0.0052 


<0.0034 


MHO 5 


1.64 


0.35 


<0.58 


<0.11* 


<0.09 


<0.23 


<0.52 


<0.45 




0.101 


0.033 


<0.034 


<0.0103* 


<0.0082 


<0.027 


<0.0083 


<0.0072 


CIDA 14 


0.27 


<0.12 


<0.12 


<0.06 


<0.06* 


<0.08 


<0.15 


<0.19 




0.040 


<0.026 


<0.019 


<0.012 


<0.013* 


<0.021 


<0.0068 


<0.0084 


KPNO 11 


<0.23 


<0.08 


<0.50 


<0.09 


<0.12 


<0.13 


<0.40 


<0.34 




<0.0073 


<0.0040 


<0.0145 


<0.0039 


<0.0056 


<0.0072 


<0.0029 


<0.0025 


V927 Tau 


<0.17 


<0.16 


<0.41 


<0.11 


<0.11 


<0.07 


<0.25 


<0.28 




<0.081 


<0.108 


<0.19 


<0.069 


<0.064 


<0.052 


<0.037 


<0.040 


MHO 7 


<0.22 


<0.07 


<0.38 


<0.20 


<0.16 


<0.20 


<0.28 


<0.13 




<0.0175 


<0.0086 


<0.032 


<0.021 


<0.018 


<0.027 


<0.0073 


<0.0033 


CFHT 7 


<0.24 


<0.18 


<1.09 


<0.15 


<0.16 


<0.27 


<0.58 


<0.31 




<0.0042 


<0.0048 


<0.0156 


<0.0036 


<0.0041 


<0.0091 


<0.0016 


<0.0009 



'^Equivalent width and line flux format as in Table 4. 

B=blend, *=marginal detection, ID=questionable line identification 



Tables 9 and 10 (continued for HST Sample) 



Star 


Ca II K 


Ca II H'^ 


[O I] A5577 


[S I] A4069 + 4076'' 


DF Tau 


10.4 


30 


0.7 


< 1.5 


Sep. 99 


31 


92 


6.6 


< 4.7 


DF Tau 


28 


38 


0.9 


3.0 


Aug. 02 


42 


56 


4.1 


4.7 


RU Lupi 


50 


39 


< 0.3 


9.2 




500 


330 


< 3 


72 


DG Tau 


61 


44 


0.39 


19' 




33 


22 


0.67 


9.2' 


TW Hya 


9.6 


32.0 


< 0.15 


< 0.8 


May 00 


65 


230 


< 2.5 


< 6.0 


TW Hya 


11.5 


45 


< 0.15 


< 0.7 


Jul. 02 


65 


270 


< 2.0 


< 4.0 


V836 Tau 


12.9 


19.5 


< 0.26* 


4.1 




1.37 


2.6 


< 0.24* 


0.79 


V819 Tau 


4.1 


4.2 


< 0.5 


< 1.7 




0.95 


1.2 


< 1.2 


< 0.8 


LkCa 7 


7.9 


7.0 


< 0.5 


< 2.5 




3.8 


4.2 


< 2.0 


< 2.0 



"^Equivalent width and line fiux format as in Table 4. 

B=blend, *=marginal detection, ID=questionable line identification 

■^Blended with He 

Lines are marginally resolved and are coadded here 
^Emission extended by 0''3 in slit and included in extraction 
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TABLE 11 
Model Accretion Parameters 



Star 


Template 


T 




L 


5 


''3600 


7*4000 




Mace 






K 


1014 cm-3 


10'^ cm % 






10-4 Lq 


10-10 Mm yr-i 






Low-mass accretors in 


Keck I sample 










vyZi lau 


10500 


0.8 


10 


1.01 


70 


A 




1 on 


tjiVL lau 


ivirivj ( 


9000 


0.7 


20 


0.53 


110 




Q1 




2M0414+2811 


CFHT 7 


7400 


2.5 


2 


3.3 


200 


20 


32 


9.0 


S0518+2327 


MHO 7 


10500 


2.5 


1 


0.59 


68 


6.7 


28 


5.0 


MHO 6 


V927 Tau 


7800 


1.5 


8 


0.31 


3.4 


0.29 


14 


4.4 


S0439+2336 


V927 Tau 


8300 


2.0 


3 


0.12 


1.4 


0.15 


6.5 


1.5 


2M0436+2351 


V927 Tau 


(9500) 


1.3 


1. 


1.7 


5.0 


0.35 


3.8 


0.53 


2M0441+2534 


CFHT 7 


8000 


1.5 


5 


0.38 


40 


3.9 


0.67 


0.19 


2M0444+2512 


CFHT 7 


8200 


2.0 


2.5 


0.034 


3.5 


0.20 


0.26 


0.092 


2M0439+2544 


CFHT 7 


(7500) 


(1.0) 


20 


0.19 


3.4 


0.14 


0.17 


0.069 


2M1207-3932(Feb.) 


CFHT 7 


7500 


0.3 


7 


1.00 


26 


0.8 


0.061 


0.023 


2M1207-3932(Nov.) 


CFHT 7 


6500 


1.5 


20 


1.4 


21 


0.7 


0.044 


0.017 


Possible accretors in Keck I sample 


2M0455+3028 


V927 Tau 


(9500) 


(2.0) 


1.5 


0.007 


0.4 


0.06 


1.7 


0.58 


MHO 5 


CFHT 7'' 


(8500) 


(2.0) 


2 


0.006 


0.22 


0.02 


0.46 


0.22 


CIDA-14 


MHO 7 


(9500) 


(2.0) 


5 


0.002 


0.14 


0.013 


0.29 


0.14 


KPNO 11 


MHO 7" 


(8500) 


(2.0) 


6 


0.008 


0.54 


0.04 


0.23 


0.062 


V927 Tau 


MHO 7 


(9500) 


(2.0) 


3 


0.001 


0.24 


0.05 


3.4'* 


1.3'' 








Accretors 


in HST sample 










DF Tau (Sep. 99) 


LkCa 7 


8500 


3.0 


2.0 


7.0 


16 


4 


2400 


460 


DF Tau (Aug. 02) 


LkCa 7 


9000 


1.6 


4.0 


2.5 


10 


1.8 


1200 


230 


RU Lupi 


V819 Tau 


9200 


2.0 


1.7 


28 


75 


13 


3500 


180 


DG Tau 


V819 Tau 


8200 


2.2 


12 


19 


12 


2.5 


870 


41 


TW Hya (May 00) 


V819 Tau 


9000 


2.0 


1.5 


12 


12 


2.2 


360 


15.4 


TW Hya (Jul. 02) 


V819 Tau 


9000 


1.8 


1.0 


12 


12 


1.7 


340 


14.6 


V836 Tau 


V819 Tau 


(9000) 


(2.0) 


0.5 


1.4 


0.65 


0.1 


32 


1.5 



"The photospheric template of V927 Tau is adjusted slightly by adding the specrum of MHO 7 

''Increased the flux of the template CFHT 7 by assuming Ay = 1 to account for the spectral type discrepancy. 

'^Decreased the flux of the template MHO 7 by assuming Ay = 0.3 to account for the spectral type discrepancy. 

''Total Lace and M for both components. 

^From BricenoetaL (2002) 



TABLE 12 

Error Estimates 



Parameter 


Uncertainty 


Dependent Variables 




Systematic 


Av 


0.5 mag 


Lace 


0.4 


7 


SpT 


0.5 classes 


M, R 


0.14 


No? 


Pre-main sequence track 


0.03-0.4 dex 


M 


-0.03 to -0.4 


Yes 


Mass calibration 


0.2 dex 


M 


to -0.2 


Yes 


Distance 


15 pc 


^acc ^ 


0.13 


No? 


Exclusion of lines 


0.25Lacc 


Lace 


+0.03 to +0.25 


Yes 


Flux Calibration 


0.18 dex 


OC LaccR^^^ 


0.26 


Sometimes 


Bolometric Correction 


0.3 dex 


Lace 


+0.3 


Yes 


Truncation Radius 


0.3 dex 




0.08 


7 


Photospheric Template 


0-0.3 dex 


Lace 


to -1-0.3 


No? 


Single/binary 


0.15 dex 


R 


-0.3 


Yes 


Re-radiation of emission 


0.3 dex 


Lace 


+0.3 


Yes 



"+ or — indicatesthat the uncertainty is only in one direction. 



TABLE 13 

Estimated M from V-band excess" 



Star 


V 


M 


R 


rv 


Av 


/5500 


Lace 


M 


Star 


mag 


Mq 


R 




mag 


10-1^ erg cm-2 s-i AA-i 


10-4 Lq 


10-10 Mq yr-1 


GM Tau 


17.577 


0.15 


0.93 


1.7 ±0.9 


2.0 


13 


79 


19 


KPNO 12 


23.228 


0.014 


0.14 


5+3 


1.6 


0.067 


0.41 


0.16 


KPNO 4 


24.722 


0.01 


0.42 


1.6 ±0.8 


4.0 


0.11 


0.67 


1.1 


KPNO 6 


22.292 


0.02 


0.68 


1 1+0.9 
1-1-0.7 


0.7 


0.044 


0.27 


0.36 



"Magnitudes, Excess, and stellar properties (except for GM Tau) from j<raus_et_^l^ (2006). 
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TABLE 14 

Comparing M measured from UV Excess and Ha 







Ha Profile" 




UV Excess'' 




Star 


Ha 10% width 


Ha EW 


logM 


Ha EW 


log M 


Ref 




(km s^"*") 


(A) 


Hog A'/f7^ vr 


(A) 


doe. Mm vr^-"-) 




iViXlW 




-ou 


1 n ft 

-±U.O 


-43 


-10.7 


1 


MHO fi 


309 


-25 


-10.3 


-55 


-9.3 


\ 


CIDA 14 


289 


-34 


-10.3 


-10 


-10.9 


\ 


ZiiVlUrt J- rtn^-^O J. -L 


362 


-250 


-10 


9Sft 
-zoo 


-o.y 


2 




215 


-102 


-11.3 




-11 1 


2 




256 


-234 


-11.3 


-600 


-10.5 


2 




320 


-24.5 




-9.5 


< —10.5 


2 


TW Hya 


361 




-9.3 




-8.7 


3 


KPNO Tau 12 


206 


-66.8 


-11.4 




-10.8-' 


2 


KPNO Tau 4 


192 


-38.4 


< -12 




-lO.O^'S 


2 


KPNO Tau 6 


235 


-41.1 


-11.4 




-11.0-f 


2 


CIDA 1 


375 


126 




232 


-7.9 


4 


GM Tau 


370 


113' 




281 


-8.7 


4 


2M1207-3932 


280 






200,400 


-11.7 


5,6 


KPNO 11 


182 


17.7'' 




13.6 


< -11.2 


5 


V927 Tau 


290 


7.0 




10.7 


< -9.9 


1 


MHO 7 


116 


9.0 




7.4 




1 


1: MuzeroUe et al. (2003) 




4: White & 


Basri (2003) 







MuzeroUe et ah 



2005) 

MuzeroUe et al. (2000) 
"Measured during Ho M measurement 
''Measured during UV excess M measurement 
'^May differ from Table 10, M, and i?* adopted 
■'^Calculated from V-band excess, see §5.1 
3M,R from Table 14 



^ohant^_et_aL 

Schaefer et al. 



(2005) 
f2006'l 



here from MuzeroUe et al. 



TABLE 15 

Line Luminosity-Accretion Luminosity Relationships for selected lines 



Diagnostic 


A (A) 


a" 


b" 






He I 


5016 


3.3 ±0.3 


1.04 ±0.05 


18 


0.30 


He I 


5876 


5.3 ±0.7 


1.46 ± 0.12'' 


12 


0.34 


He II 


4686 


3.7 ±0.5 


1.04 ±0.08 


16 


0.38 


Nal D 


5893 


3.3 ±0.7 


1.09 ± 0.11 


11 


0.41 


Call 


8542 


2.5 ±0.5 


1.02 ± 0.11 


20 


0.44 


OI 


6300 


2.8 ±0.8 


0.96 ± 0.16 


26 


0.67 


Ha 


6563 


2.0 ±0.4 


1.20 ±0.11 


35 


0.57 


H/3 


4861 


2.6 ±0.2 


1.22 ± 0.05 


19 


0.26 


H7 


4340 


3.0 ±0.2 


1.24 ±0.04 


19 


0.21 


Call K 


3932 


2.0 ±0.3 


0.96 ± 0.05 


19 


0.33 


H-9 


3771 


3.4 ±0.2 


1.17 ±0.03 


19 


0.18 


" log iiine = 


a + 61c 


giacc, with Lace and L;i„e in 


units of Lq 



''Number of accretors with detections. 
'^Standard deviation of detected line fluxes from linear fit. 
''Significant deviation from 1.0, use with caution outside our M range. 
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Wavelength (A) 



5500 



Fig. 1. — Three spectra at 3300-5600 A. The spectrum of V927 Tau is dominated by photospheric emission, with strong TiO bands and 
Ca I A4227 absorption line. Moderate accretors, such as S0439+2336, are dominated by photospheric emission longward of the Balmer 
jump at 3700 A but excess Balmer continuum emission at shorter wavelengths. The TiO bands at < 5600 A axe only weakly detected from 
the strong accretor S0518+2327, which shows a heavily veiled spectrum both shortward and longward of the Balmer jump. 



Accretion onto very low-mass stars and brown dwarfs 27 




Fig. 2. — The Ca I A4227 absorption line in our LRIS spectra (solid lines) compared with a photospheric template (dotted lines, template 
listed in Table 5). When present, accretion fills in some of the absorption. The veiling is calculated from the depth of the absorption line. 
Emission lines in the vicinity of the Ca I A4227 line from several strong accretors (e.g., CIDA 1, GM Tau) limit our ability to measure high 
veilings and may lead to overestimated small veilings. 
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Fig. 3. — A TiO absorption band at 4955 A used for veiling measurements. 
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1.5 
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Fig. 4. — A TiO absorption band at 6150 A used for veiling measurements. 
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3500 3600 3700 3800 3900 4000 3500 3600 3700 3800 3900 4000 

Wavelength (A) Wavelength (A) 



Fig. 5. — The spectral region from 3500-4000 A for late-K/early M stars (upper left), M5 stars (upper right), M5.5-M6 stars (lower left), 
and M6.5-M8 stars (lower right). The Balmer continuum, Balmer series, and Ca II H & K are apparent in emission. 
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3400 3800 4200 3400 3800 4200 



Fig. 6. — The Balmcr jump region for stars with weak an observed spectrum to a photospheric 

template (dotted line, both smoothed by 11 pixels) increases at ^ 3700 when aScretion is present. The solid line shows the same ratio, after 
subtracting excess Balmcr line and continuum emission from the spectrum. For a perfect model and photospheric template, the residuals 
should be a line equal to the flux ratio with noise introduced from both the spectrum and photospheric template. Some excess Balmer 
continuum emission is seen from each source plotted here. The excess Balmer continuum emission from S0439+2336 is strong enough to 
be attributed to accretion. Fainter Balmer excess emission may be produced by either chromospheric activity or accretion. 
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Fig. 7. — Contours of acceptable T and rie for isothermal hydrogen slab model fits to the spectra of 2M0414+2811 (left) and 2M1207- 
3932 (right) . The vertical contours that determine T arc constrained by the slope of the Balmer continuum. The horizontal contours that 
determine Ue are constrained by the size of the Balmcr jump. For 2M0414+2811 the acceptable values of T and Tie together differ from 
their intersection because they are only acceptable for different values of L. 
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Fig. 8. — The observed emission (soUd spectrum), corrected for reddening (Table 1), is fit with the combination of a photospheric 
template (dashed line, listed in Table 5 and not shown in the smaller panel) and a synthetic accretion spectrum from a slab model (fight 
solid line, Balmer lines are not shown for clarity). The total fit is shown as a dashed line that overlaps the observed spectrum. 
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-13 -12 -11 -10 -9 -8 -12 -11 -10 -9 -8 -7 

log (Mq yr"') from UV Excess log (Mg yr"') 

Fig. 9. — UV excess measure of M versus non-simultaneous measurements of M from models of the Ha line profiles (left, from Muzerolle 
et al. 2000, 2003, 2005) and versus non-simultaneous measurements of the Ha: 10% width (left). The circles are accretors with M measured 
from UV excess and the diamonds are accretors with M measured from the V-band excess (Table 14) reported by Kraus et al. (2006). The 
dashed line on the left panel denotes where two accretion rate measures would be equal. The dashed line on the right panel denotes the 
M-Ha 10% width relationship calculated by Natta et al. (2004). 
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Fig. 10. The relationship between selected Hne luminosities and the aceretion eontinuum luminosity. The fiUed cireles represent the 
accretors with strong excess Balmer continuum emission that is produced by accretion. The upper hmits in Lace are the five possible 
accretors, which have weak Balmer continuum emission that may be attributed to either accretion or chromospheric activity. We include 
results from Valenti et al. (1993, empty diamonds and corrected for the 140 pc distance to Taurus) and GuUbring et al. (1998, filled squares) 
for fits to the H/3, H7, and Ca II K line luminosities, and Hartigan et al. (2003, filled triangles) for fits to the Ha, Ca II A8542, and [O I] 
line luminosities. The dashed lines are linear fits to the data as described in §5.3 and Table 15. 
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Fig. 11. — Ca II A8542 (upper left) and Ca II A8662 (upper right) line flux versus mass accretion rate. Wc compare our data (filled 
circles for accretors, empty circles for the five possible accretors) with that from Hartigan et al. (2003, filled diamonds for Ca II A8542), 
MuzeroUe et al. (1998, filled triangles for Ca II A8542 and He I A5876) and Mohanty et al. (2005, upside-down triangles for Ca II A8662). 
The empty upside-down triangles on the upper right panel signify that the M was measured by Ha modelling and is multiplied by 5. The 
dashed lines in both panels show the best fit through the data. 
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Fig. 12. — The slope of the Balmer continuum versus stellar mass (left) and M (right). The circles are from our Keck //LRIS data, the 
squares from HST/STIS G430L spectra. The slope of the Balmer continuum increases to smaller M and M, suggesting that it is produced 
at cooler gas. 
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Fig. 13. — The Balmer jump versus stellar mass (left) and A/ (right), as in Fig. 12. The squares are from iGuUbring et al.l ()1998l ). The 
Balmer jump tends to increase to smaller M and M, suggesting that the accretion continuum is produced in gas with a lower electron 
density. 
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Fig. 15. — The relationship between mass and Lace (left) and 



(right), with circles from our data and open triangles from 



IGuIIbring etaLl HlMB ). The sensitivity of and y^!^ from the UV excess improves to lower masses. 




Fig. 14. — Mass versus mass accretion rate for Taurus members with Ay < 2.5 and measured accretion rates. The data plotted are M 
measured from excess UV continuum emission in this work (filled circles), from excess UV continuum emission measured by GuUbring et 
al. (1998, open triangles), optical veiling measured by Hartigan & Kenyon (2003, asterisks), UV excess measured by Calvet et al. (2004, 
open diamonds), and Ha profile modelling (open squares, multiplied by 5, MuzeroUe et al. 2005). The dashed line shows M oc — 7.7M^ *^. 



